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ABSTRACT 

Aims. We study three Galactic H n regions - RCW 79, RCW 82 and RCW 120 - where triggered star formation is taking place. Two 
stellar population are observed: the ionizing stars of each H n region and young stellar objects on their borders. Our goal is to show 
that they represent two distinct populations, as expected from successive star forming events. 

Methods. We use near-infrared integral field spectroscopy obtained with SINFONI on the VLT to make a spectral classification. We 
derive the stellar and wind properties of the ionizing stars using atmosphere models computed with the code CMFGEN. The young 
stellar objects are classified according to their A*-band spectra. In combination with published near and mid infrared photometry, we 
constrain their nature. Linemaps are constructed to study the geometry of their close environment. 

Results. We identify the ionizing stars of each region. RCW 79 is dominated by a cluster of a dozen O stars, identified for the first 
time by our observations. RCW 82 and RCW 120 are ionized by two and one O star, respectively. All ionizing stars are early to late 
O stars, close to the main sequence. The cluster ionizing RCW 79 formed 2.3±0.5 Myr ago. Similar ages are estimated, albeit with a 
larger uncertainty, for the ionizing stars of the other two regions. The total mass loss rate and ionizing flux is derived for each regions. 
In RCW 79, where the richest cluster of ionizing stars is found, the mechanical wind luminosity represents only 0.1% of the ionizing 
luminosity, questioning the influence of stellar winds on the dynamics of these three H n regions. The young stellar objects show four 
main types of spectral features: H 2 emission, Bry emission, CO bandheads emission and CO bandheads absorption. These features 
are typical of young stellar objects surrounded by disks and/or envelopes, confirming that star formation is taking place on the borders 
of the three H n regions. The radial velocities of most YSOs are consistent with that of the ionized gas, firmly establishing that their 
association with the H n regions. Exceptions are found in RCW 120 where differences up to 50 km s~ ] are observed. Outflows are 
detected in a few YSOs. All YSOs have moderate to strong near-IR excess. In the [24] versus A"- [24] diagram, the majority of the 
sources dominated by H2 emission lines stand out as redder and brighter than the rest of the YSOs. The quantitative analysis of their 
spectra indicates that for most of them the IT emission is essentially thermal and likely produced by shocks. We tentatively propose 
that they represent an earlier phase of evolution compared to sources dominated by Bry and CO bandheads. We suggest that they still 
possess a dense envelope in which jets or winds create shocks. The other YSOs have partly lost their envelopes and show signatures of 
accretion disks. Overall, the YSOs show distinct spectroscopic signatures compared to the ionizing sources, confirming the presence 
of two stellar populations. 

Key words. ISM: H 11 regions - ISM: bubbles - Stars: formation - Stars: early-type - Stars: fundamental parameters - Stars: winds, 
outflows 

1. Introduction A particular case concerns star formation on the borders 

of H n regions. Accord ing to the collect and collapse model 
Massive stars play a significant role in several fields of as- jElmegreen & Ladalll977h . a dense shell of material is trapped 
trophysics. They produce the majority of heavy elements and between the shock and ionization front of an expanding H n re- 
spread them in the interstellar medium, taking an active part g i ons . when the amount of collected material is large enough, 
in the chemical evolution of galaxies. But they also end their global shell fragmentation occurs and new stars are formed. The 
life as supernovae and gamma-ray bursts. Through their strong observation of molecular condensations on the borders of several 
winds and ionizing fluxes they power H 11 regions and bubbles h ii regions and the subsequent identification of YSOs within 
which are often used to trace metallicity gradients in galaxies, these clumps dDeharveng et alj200 3. 2005; Zava gno etalT2 006. 



The energy they release in the interstellar medium is thought 2007J iDeharveng et al.ll2008l 120091; iPomares et all 120091) con 



to trigger second-generation star formation events. Observations fi rms ma t this mechanism is at work at least in some H 11 regions. 

of young stellar objects (YSO) in molecular clouds surrounding , , ,. , „ 

. / _ . , , ,. ,. . Other mechanisms leading to triggered star formation ex- 

(clusters of) massive stars lend support to this mechanism (e.g. . . „ . ... ,. , ° ., 7, * j n a 1 

0.. — 1 rr; m ?Z" 1, 1st. Some work qualitatively as the collect and collapse model 

Walborn et al. 2002; Hatano et al. 2006). . . c , , ■ T \ 
in the sense that the clumps are formed during the H 11 region 

expansion . For instance, dynamical instabilities of the ioniza- 



Send offprint requests to: F. Martins tion front dVishniadll983l;lGarcia-Segura & Francoll 19961) create 

* Based on observations collected at the ESO Very Large Telescope molecular condensations separated by zones of lower densities, 
(program 081.C-0057) The newly formed clumps grow until they become Jeans unsta- 
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ble and collapse. Alternatively, second generation star formation 
can happen in pre-existing clumps. If the neutral gas in which 
the H n region expands is not homogeneous, the outer layers 
of the molecular overdensities are ionized like the borders of a 
classical H n region. A shock front precedes the ionization front 
inside these clumps, lead ing to their collapse dDuvert et al J 1990t 
iLefloch & Lazarefflll994l) . 

A number of questions regarding triggered star formation 
remain unanswered. The properties of the observed YSOs are 
poorly known besides a crude classification in class I or class 
II objects by analogy with low-mass stars. YSOs usually dis- 
play near-infrared s pectra with CO, Bry and/or H2 emission 
lines dBik et al1l2006h . The relation, if any, between objects with 
different spectroscopic appearance is not clear. Besides, in the 
regions where the collect and collapse process is at work, the 
quantitative properties of the ionizing sources of the H 11 re- 
gions are not known. In particular, the relative role of ionizing 
radiation and stellar winds on the dynamics of such regions is 
debated. The timescales under which material accumulates and 
fragments, and the properties of t he resulting clumps d epend on 
the strength of those two factors dWhitworth et aljl994t) . Hence, 
one might wonder whether the nature of the newly formed ob- 
jects depends on the properties of the stars powering the H n 
regions. 

In the present study, we tackle these questions by investi- 
gating the properties of the ionizing stars and YSOs of three 
Galactic H n region s: RCW 79, RCW 82 and RCW 120 
(Rodgers et al.l Il960h . Those regi ons are known to be the 
sites of triggered star formation dZavagno et al] 120061 [2007; 
iPomares et al]l2^lKanglt^l2009l) . We have used SINFONI 
on the VLT to obtain near-infrared spectra of both the ionizing 
stars and a selection of YSOs in each region. Our main goals 



- Identify the ionizing stars of each region, and derive their 
stellar properties using atmosphere models. In particular, we 
want to determine their ionizing fluxes and mass loss rates 
in order to better understand the dynamics of the H 11 re- 
gions. Equally important is the determination of the age of 
those stars, since it can be related to the presence of YSOs to 
quantitatively confirm the existence of triggered star forma- 
tion. 

- Constrain the nature of the YSOs on the borders of the 
H 11 regions. In combination with infrared photometry, spec- 
troscopy can reveal the presence of disks or envelopes. The 
evolutionary status of those objects can thus be better under- 
stood. In particular, it can be clearly shown whether they are 
stars still in their formation process or objects already on the 
main sequence. 

The paper is organized as follows. In Sect. [2] we present the 
three H 11 regions targeted in this study. Sect. [3] describes our 
observations. The analysis of the ionizing stars is presented in 
Sect.|H while the YSOs are discussed in Sect. [5] We discuss our 
results in Sect. [6] and summarize our conclusions in Sect. [7] 



2. Presentation of the observed H n regions 

2.1. RCW 79 

RC W 79 is a Southern Galactic H 11 region lo cated at 4.2+1 
kpc (Russeil et al. 1998). Its diameter is ~6.4 pc. Zavagno et al. 
(2006j) (hereafter ZA06) used Spitzer GLIMPSE and SEST- 
SIMBA 1 .2-mm continuum data to study the star formation on 



the borders of this region. A layer of warm dust is clearly seen 
at 8 fim surrounding RCW 79 (Fig. [TJ. A compact H n region 
is observed in this layer as well as five cold dust condensa- 
tions (masses between 100 and 1000 M ) detected by 1.2-mm 
continuum emission. YSOs have also been revealed by Spitzer- 
GLIMPSE observations, leading ZA06 to conclude that trig- 
gered massive star formation was at work possibly thro ugh the 
collect and collapse process (Elmegr een & Ladal Il977l) . Eight 
YSOs detected in the near IR have been observed with SINFONI 
at the VLT. They are identified on Fig.Q] In several cases, mul- 
tiple sources have been uncovered by the SINFONI observa- 
tions reported in this paper. They are identified in Sect. [5] and 
Appendix iBl 

The ionizing stars of RCW 79 were unknown before the 
present study. Using photometry and choosing objects inside the 
H n region contours, we have identified and observed a cluster 
of possible OB stars candidates (blue square in Fig . [TJ . They are 
described in Sect. [4] 




-bl:U:l» 



Fig. 1. Identification of the RCW 79 YSOs targeted in the 
present study on a composite infrared color image. Blue and 
green are the Siptzer-GLIMPSE images at 3.6 and 8.0 //m re- 
spectively. Red is an unsharped masked image obtained from 
the Spitzer-MIPS frame at 24 /mi (as explained in ZA07). The 
bright extended emission has been subtracted to emphasize the 
YSOs. The candidate ionizing sources are located inside the blue 
rectangle (see also Fig.|4]i. 



2.2. RCW 82 

RCW 82 is a Southern H n reg ion located at a distanc e of 3.4+1 
kpc. Its radius is about 5 pc. Pomar eTet al.l d2009) (hereafter 
PO09) detected molecular material in 12 CO(1-0) and 13 CO(1-0) 
surrounding the H 11 region. They showed that some of the struc- 
tures correspond to dense material collected between the shock 
front and the ionization front during the expansion of the H 11 re- 
gion. Masses of these clumps range from 200 to 2500 M . Star 
formation is observed on the borders of RCW 82, with a total of 
63 candidate YSOs. Among these, we have selected five YSO 
candidates visible in the A'-band for our study. These objects are 
shown on Fig. [2] 

PO09 identified four candidate ionizing sources, among 
which two are likely O stars (see their Fig. 5). These sources 
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are located just south of the 24 fim emission ridge in Fig [2] and 
better displayed in the upper panel of Fig. [5] 
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Fig. 2. Same as Fig.Q]for RCW 82. 



2.3. RCW 120 

A deta iled study of RC W 120 can be found in Zavagno et al.l 
d2007l) (here after ZA 07). iDeharveng & Zavagnol d2008b (DE08) 
and lDeharvengetalJ d2009j) (DE09). RCW 120 is the nearest of 
the three H n regions observed in the present study. Its distance, 
D= 1.35+0. 3 kpc, is well determined since both the photometric 
and kinematic distances are in good agreement (DE09). 

RCW 120 has a circular geometry (diameter -3.5 pc) and 
is surrounded by a shell of dense material collected during the 
expansion of the H n region. The cold dust emission of this shell 
has been observed with the ESO SEST at 1.2-mm (ZA07) and 
with the APEX-LABOCA camera at 870 //m (DE09). Its mass 
is in the range 1100-1900 M ; it is fragmented with massive 
fragments elongated along the ionization front. Star formation is 
at work in these condensations, as discussed by ZA07 and DE09. 
Twelve candidate YSOs have been observed with SINFONI and 
are identified in Fig. [3] 

The central ionizing star of RCW 120 is CD-38°11636. 
Its spectral type (estimated from spectrograms) is 08 
dGeorgelin & Georgelinl 1 1970b or 09 dCramptonl Il97ll) - see 
also the discussion in DE08. I ts extinction, determined by 
lAvedisova & Kondratenkol d 19841) is Ay =4. 65 mag. It is located 
just north of the bright 24 //m emission ridge in Fig. [3] Its char- 
acteristics are discussed in Sect. [4] 



3. Observations and data reduction 

For the current project, we selected two types of sources: the 
candidate ionizing stars of the H n regions and a few candidate 
YSOs. For the latter, we selected the brightest objects in the K- 
band. 

Data have been collected at the ESO-VLT on April 25* 
and 26 th 2 008. The near infrared integral field spectrograph 
SINFONI dEisenhaueretalJl2003l) was used in seeing-limited 
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vV* : 28:00 
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Fig. 3. Same as Fig. Q]for RCW 120. On this image, the ex- 
tended high brightness emission has been subtracted to enhance 
the small scale structures and the point sources. 



mode to obtain medium resolution /f-band spectra of our se- 
lected sources. We selected the 250 mas scale which provided 
a field of view of 8"x8". Sequences of source and sky expo- 
sures were conducted in order to ensure optimal sky subtraction. 
For the faintest sources, two exposures were made, with a 1" 
offset between them to minimize pixels artifacts. The exposure 
times ranged from 1 minute for K-9 sources to 10 minutes for 
K-13 objects. Telluric stars were observed regularly during the 
night to allow proper atmospheric correction. Standard calibra- 
tion data were obtained by the ESO staff. The observations were 
conducted under an optical seeing ranging from 0.6" to 1.2". 

Data reduction was made with the SPRED software 
dAbuter etalj 12006). After bias subtraction, flat field and bad 
pixel corrections, wavelength calibration was done using a Ne- 
Ar spectrum. Fine tuning was subsequently performed using sky 
lines. Telluric correction was done using standard stars from 
which Bry and, when present, Hei 2.112 /im were removed. 
The resulting spectra have a signal to noise of 10-100 depend- 
ing on the source brightness and wavelength. Their resolution is 
~ 4000 (see also Sect. 15.41 ). They were extracted with the soft- 
ware QFitsView Q, carefully selecting individual "source" pix- 
els one by one to avoid contamination by neighboring objects. 
Background pixels were selected close to the "source" pixels in 
order to correct for the underlying nebular emission. In practice, 
the average spectrum of these background pixels is subtracted 
to the each individual source pixels. Then all corrected source 
pixels are added together to ensure optimum extraction of the 
pure stellar sprectrum. This method is possible with integral field 
spectroscopy and was one of the main drivers for the choice of 
SINFONI for the present observations. 



1 http ://w w w. mpe. mpg . de/~ ott/QFi ts View/ 
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Fig. 4. A'-band images of the ionizing star cluster of RCW79 
from 2MASS (left) and SINFONI (right). The intensity of the 
SINFONI image corresponds to the median of the spectral di- 
mension of the datacube. 



4. Ionizing stars 

In this section we describe the spectral morphology of the ioniz- 
ing sources of each region. We determine their stellar and wind 
properties using atmosphere mod els computed with the code 
CMFGEN dHillier & MilleJll998l) . Fig. shows the cluster of 
stars responsible for the ionization of RCW 79. The stars are la- 
beled on the 2MASS /f-band image presented on the left part 
of the figure. Numbers correspond to decreasing 2MASS K— 
band magnitudes. Note that star number 3 is not indicated: it 
corresponds to a source out of the field which turned out to be 
a foreground star (its spectrum displays CO bandheads in ab- 
sorption typical of cool, evolved stars). The right panel of Fig. [4] 
shows a mosaic of the SINFONI fields observed. In many cases, 
the higher spatial resolution reveals several components for the 
same 2MASS source (e.g. source 9). Consequently, we have as- 
signed new names to the resolved components. The identifica- 
tion is shown on the right side of Fig. |4] Similarly, the ionizing 
sources of RCW 82 and RCW 120 are displayed in Fig. 



4.1. Photometry and spectral classification 

Table Q] summarizes the observational properties of the ionizing 
stars of the three regions. Photometry is from 2MASS. When a 
given 2MASS point source was resolved in several components 
by our observations, the relative SINFONI fluxes were used to 
recompute the individual magnitudes, using the 2MASS magni- 
tude as a measure of the total flux from all components. In ad- 
dition, the /f-band magnitude of star 2 was revised because its 
magnitude is clearly different from that of star 1, in contradiction 
to the the 2MASS values (see right panel of Fig.|4j>. Inspection of 
the 2MASS catalog revealed that the aperture used to compute 
the /f-band magnitude was wider than the separation of sources 
1 and 2. Hence, the flux used to estimate the /f-band magni- 
tude of star 2 was contaminated by flux from star 1 . We used 
the SINFONI flux and 2MASS magnitude of star 1 to calibrate 
the /f-band magnitude of star 2 from the observed SINFONI 
flux. The //-band magnitudes are from 2MASS. Since there is 
no SINFONI //-band data, one cannot recalculate the magnitude 
of the components of unresolved 2MASS sources. In that case, 
we give the 2MASS //-band magnitude as a lower limit. 
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Fig. 5. Same as Fig.fflfor RCW 82 (top) and RCW 120 (bottom). 
North is up and East is left. 



For each region, we calculated the absolute Zf-band magni- 
tude using information on the distance and extinction: 

- RCW 79: a dista nce of 4.2+1 kpc was derived by 
iRusseil et al.l d 1 998b . We estimated the Zf-band extinction 
from the H - K color of the early type stars of the region 
(see below for the spectral classification). Since (H - K)o 
is b asically independent of the spectral type for O stars 
( see iMartins & Plezll2006l) . one can use the 2MASS magni- 
tudes even in the cases were SINFONI revealed several sub- 
components. In practice, we used stars 4, 5, 7, 8, 9 and 10 to 
obtain A K =0.7 ±0.1. 

- RCW 82: We proceeded as for RCW 79 to calculat e MK as- 
suming a distance of 3.4+1 kpc (IRusseil et al.l 1998b . The ex- 
tinction estimate was based on stars RCW82 e2 and RCW82 
e3. Both of them are O stars. A mean value of Ak=0.36±0.2 
was derived. 

- RCW 120: With a distance of 1.35+0.3 kpc (ZA07), 
RCW 120 is the closest region of our sample. The region 
is ionized by a single O star. The estimated extinction is 
Ak=0.50 ±0.1. This is in good agreement with previous es- 
timates (see Sect. 12.31 ). 



Figs I A. 1 1 and I A. 21 show the spectra of the candidate ionizing 
stars of RCW 79, RCW 82 and RCW 120. Most of the stars show 
Bry and He 1 2. 1 12 yi/m absorption and, depending on the sources, 
Hen 2.189 fim absorption, Civ 2.070-2.083 fj.m emission and 
the N m/C m/O in emission complex at 2. 1 15 yt/m. These features 
are ty pical of OB stars. We relied on the atlases of Zf-band spec- 
tra of Hanson et al. (1 19961 120051) to assign spectral types to the 
ionizing stars of RCW 79, RCW 82 and RCW 120. In practice, 
we have used the following scheme: 

• the presence of He n 2. 189 yum indicates a spectral type ear- 
lier than 08 

• C iv 2.070-2.083 fim emission is observed in 04-6 stars and 
is the strongest at 05 
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• the N m-C m-O m emission complex is observed in emission 
together with He i 2.1 12 //m in absorption in 04-07 stars 

• He i 2. 1 12 yum disappears at spectral types later than B2.5 

• 03-7 supergiants have Bry either in emission or in absorp- 
tion weaker than He n 2. 1 89 /urn. Later type supergiants have 
a narrow Bry absorption profile from which Hei 2.164 yimi 
is clearly separated. Dwarfs and giants have broader Bry ab- 
sorption (He i 2. 164 yt/m is blended with Bry). 

Based on these features, we classified the ionizing stars as OB 
dwarfs or giants. The results are listed in Table Q] The earli- 
est spectral types are 04-6 and are observed only in the rich- 
est ionizing cluster (RCW 79). Most stars are late O - early B 
dwarfs/giants. We do not observe any supergiant. The distinction 
between dwarfs and giants is not easily feasible with our spectra. 
Hence we adopt the conservative approach to give the luminos- 
ity class V/III to most of our stars. Some of the stars did not 
show any of the above features. Instead, strong CO bandheads 
were clearly visible. Such spectra are typical of cool, evolved 
low mass stars Q. They do not contribute any ionizing flux and 
are most likely foreground or background stars. We classify them 
as "late" in Table Q] Three stars also show strong Hei 2.056 
fim and Bry emission, without Hei 2.112 yimi emission. In one 
case (RCW82 el), the Hi Pfund serie and Mgn 2.138-2.144 
/ um we re also de tected in emission. A ccording to lHanson et alj 
(1996) (see also IClark & Steele! l2000j). these features are seen 
in Oe-Be stars. Such objects are thought to host a circumstel- 
lar disk contributing a significant fraction of the near-IR emis- 
sion (both line and continuum). Note that the ionizing source of 
RCW 120 seems to be double (Fig. [5]). We tried to extract the 
spectra of both components, but found no difference between 
them. This indicates either that both components have the same 
spectral type, or that the spatial resolution of our observations is 
not sufficient to separate the spectra. Hence, we have treated the 
ionizing source of RCW 120 as a single source. 

PO09 estimated a spectral type 06.5V and 07.5V for stars 
RCW82 e2 and RCW82 e3 respectively from J vs J - K and 
/ - H vs H - K diagrams. Our spectroscopic classification in- 
dicates later types: 09-B2V/III for both stars. Given the diffi- 
culty to make a spectral classification from near-infrared color- 
color and color-magnitude diagrams, the difference is accept- 
able. Similarly, star RCW82 e4 is identified as a late-type giant 
which is confirmed by our spectroscopy. PO09 showed that star 
RCW82 el presented a near-IR excess, which is consistent with 
the presence of a disk in a star of spectral type Oe or Be. ZA07 
reported a spectral type 08V for the ionizing star of RCW 120. 
We prefer a spectral type 06-8 V/III, in rather good agreement. 
Finally, the ionizing sources of RCW 79 were not previously 
identified. ZA06 estimated that a single 04V star could power 
the H ii region. This is roughly consistent with our finding that a 
cluster of about ten 04 to 09 stars is responsible for the ioniza- 
tion. 

4.2. Spectroscopic analysis 

In this section we derive the stellar and wind properties of the 
ionizing stars of the H n regions by means of spectroscopic anal- 
ysis with atmosphere models. We subsequently use the derived 
properties to constrain the age of the ionizing populations. 



2 CO absorption is also observed in red supergiants. However, those 
stars are much brighter than O stars in the A"-band. This is not the case 
of the stars with CO absorption of our sample. 



4.2.1 . Stellar and wind properties 

The stellar and wind properties have been derived through spec- 
troscopic analys is. Atmosphere models were computed with the 
code CMFGEN dffiUier & Miller]l998l) . This code solves the ra- 
diative transfer and statistical equations in the co-moving frame 
of the expanding atmosphere, for light elements as well as for 
metals. It thus produces non-LTE, line blanketed atmosphere 
models with winds. An exhau stive description of th e code and its 
approximations is given in Hillie r & Millerl dl998l). Summaries 
of the main characteristics can also be found in iHillier et al.l 
d2003h : lMartins et alj d2008l) . The resulting synthetic spectra are 
compared to the SINFONI /^-band spectra to constrain the main 
physical parameters. In practice, we have proceeded as follows: 

• Effective temperature: the determination of T e ff usually relies 
on the ratio of He i to He n lines. We have used the following 
lines: Hei 2.112 /mi and Hen 2.189 /mi when present. Crv 
2.07-2.08 yum lines are used as secondary in dicators since 
they a ppear in emission in the hottest O stars dHanson et al.l 
1996). An uncertainty of 2000 to 3000K is achieved depend- 
ing on the star. For stars cooler than ~32000 K Hen 2.189 
//m vanishes. The uncertainty on T e ff is large (5000 K) and 
is set by the presence and strength of He i lines. 

• Mass loss rate: The main indicator in the /f-band spectrum 
of O stars is Bry. It is filled by emission as the wind strength 
increases. We have used this line to constrain the mass loss 
rates of the ionizing stars. However, as demonstrated below, 
it becomes almost insensitive to M below ~ 3 x lO _8 M yr~ 1 
so that we could derive only upper limits on M. We stress 
once more that our correction for nebular contamination en- 
sures the best extraction possible of the stellar line profiles. 
Thus, our mass loss rate determination does not suffer from 
major uncertainties due to non stellar emission. 

• Luminosity: log y- was derived from the absolute magnitude 
and the K-h&nA bolometric correcti on. The latter was com - 
puted from T e g and the calibration of Martin s~& Plezl ((2006). 
The uncertainty on log y- was calculated from a full error 
propagation and takes into account the uncertainties on the 
distance, T e g and A K . It is of the order 0.2 dex for stars in 
the three regions. 

Several parameters could not be derived due to the lack 
of spectral diagnostics. In particular, the surface gravity was 
adopted. Since most of our stars are dwarfs/gia nts, we chose 
logg = 4.0 which is typical of O dwarfs dMartins et al.l 
l2005al). In principle Bry can be used to constrain logg (e.g. 
iRepolust et al.ll2005l) . However, the S/N of a few tens combined 
to the medium resolution of our spectra prevents an accurate de- 
termination. Besides, since we use this line to determine mass 
loss rates and since the line core depends not only on M but also 
on logg (although in a less sensitive way), we decided to fix 

Modern atmosphere models for massive stars also in- 
clud e clumping si n ce evidence for wind inhomogeneities exis t 
(e.g. IHillier! [19911 IHillier etalJ 12001 iLeoine & Moffatl fl999l) . 
Clumping is usually quantified by a volume filling factor /. 
Currently, no diagnostics of clumping have been identified in 
the K-band spectra of O stars with moderate winds. Such di- 
agnost ics are traditionally found in the UV, optical and submm 
range dHillierill99ll:lEversberg et alJll998l:lBlomme et al.ll2002l) . 
Hence, we have simply adopted the canonical v alue of 0. 1 for / 
dHamann & Koesterke1ll998t IHillier et aill2001l) . 

Another important parameter is the wind terminal velocity 
(v ra ). It is determined from the blueward extension of UV P- 
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Table 1. Position and photometry of ionizing stars. 



Source 


RA 
h m s 


DEC 

°' " 


H 


K 


MK 


ST 


RCW79 1 


13:40:09.11 


-61:43:51.9 


9.10+0.04 


8.29±0.03 


-5.53±0.53 


OBe 


RCW79 2 


13:40:09.64 


-61:43:50.1 


> 10.59 


10.41±0.06 


-3.46±0.53 


07.5-8V/III 


RCW79 4a 


13:40:12.63 


-61:44:16.0 


>9.41 


9.23±0.02 


-4.64±0.53 


04-6V/III 


RCW79 4b 


13:40:13.10 


-61:44:17.1 


>9.41 


11.29±0.02 


-2.54±0.53 


? 


RCW79 5a 


13:40:12.19 


-61:43:47.6 


>9.40 


9.24±0.05 


-4.62±0.53 


04-6V/III 


RCW79 6a 


13:40:11.47 


-61:43:30.6 


>9.24 


9.65±0.07 


-4.22±0.53 


O6-8V/III 


RCW79 6b 


13:40:10.98 


-61:43:29.6 


>9.24 


11.65±0.07 


-2.22±0.53 


09-B2V 


RCW79 6c 


13:40:11.17 


-61:43:30.9 


>9.24 


12.97±0.07 


-0.86±0.53 


>B2.5V 


RCW79 7 


13:40:12.93 


-61:43:47.8 


9.93±0.04 


9.64±0.03 


-4.23±0.53 


O6-8V/III 


RCW79 8a 


13:40:11.98 


-61:43:32.9 


>9.93 


10.33±0.04 


-3.54±0.53 


09-B2V/III 


RCW79 8b 


13:40:11.82 


-61:43:32.7 


>9.93 


10.55±0.04 


-3.32±0.53 


09-B2V/III 


RCW79 9a 


13:40:13.06 


-61:43:30.6 


>10.12 


10.79±0.03 


-3.07±0.53 


09-B2V/III 


RCW79 9b 


13:40:12.95 


-61:43:28.6 


>10.12 


10.91±0.03 


-2.96±0.53 


09-B2V/III 


RCW79 9c 


13:40:12.56 


-61:43:29.5 


>10.12 


11.87±0.03 


-1.95+0.53 


B 


RCW79 10 


13:40:12.56 


-61:43:58.9 


10.09±0.07 


9.84±0.05 


-4.05±0.53 


06.5-8V/III 


RTW7Q 1 Onp 

JXV^ VV / 7 1L/I1C 


1 3 -4.0 1 3 44 


KJl .tJ.Jj.7 


19 71+0 OQ 






Be 


RCW79 lOwl 


13:40:12.06 


-61:43:57.4 


>11.69 


12.11+0.05 


-1.71+0.53 


B0-2.5V 


RCW79 10w2 


13:40:11.79 


-61:43:57.8 


>11.69 


12.04±0.05 


-1.78+0.53 


B0-2.5V 


RCW79 10w3 


13:40:11.19 


-61:44:00.1 


12.31+0.10 


11.99±0.07 


-1.83±0.53 


Bl-2 


RCW79 1 1 


13:40:13.44 


-61:43:41.0 


10.39±0.04 


9.93±0.03 


- 


late 


RTW7Q 1 ?a 

IVV^ VV 17 l^d 


1 "V40- 1 ? 4? 


-(S1 -44-04 ? 


> 10.67 


1 n is+o 04 




late 


RCW79 12b 


13:40:12.19 


-61:44:04.9 


> 10.67 


11.49±0.04 


-2.38±0.53 


09-B2V/III 


RCW79 CHII 6 


13:40:53.71 


-61:45:46.8 


11.55±0.07 


10.72±0.05 


-4.03±0.53 


07.5-9.5V/III 


RCW82 el 


13:59:29.86 


-61:23:08.9 


9.55+0.03 


9.09±0.02 




Be 


RCW82 e2 


13:59:29.13 


-61:23:04.2 


9.59±0.08 


9.40+0.03 


-3.62±0.67 


09-B2V/III 


RCW82 e3 


13:59:28.09 


-61:23:00.4 


9.41±0.03 


9.19±0.03 


-3.83±0.67 


09-B2V/III 


RCW82 e4 


13:59:27.37 


-61:23:20.3 


10.44±0.02 


9.35±0.02 




late 


RCW120e 


17:12:20.82 


-38:29:25.5 


7.71±0.04 


7.52+0.02 


-3.62±0.49 


O6-8V/III 



in the fitting process. Since it is proportional to M/v^ an error on 
Voo translates into an error on M. However, it is very unlikely that 
our estimates of v'co are off by more than a factor of 2 (which is 
already a rather extreme case). Hence, a rather conservative up- 
per limit of 10~ 7 M yr (including the uncertainty on v ra and 
on the normalization and S/N of the observed spectrum) is deter- 
mined. Since star RCW79 5a is one of the earliest of our sample, 
it is expected to show the strongest mass loss rate. Consequently, 
the upper limit on M derived for this star also applies to all other 
ionizing stars. The upper limits we derive are fully consistent 
with recent determination of mass loss rate s of Galactic Q5-9 
dwarfs based on opti cal and UV diagnostics dBouret et al.| [2Q03: 
Repo lust et alj 120041) . Actually, mass loss rates up to two or- 
ders of magnitudes smaller are routinely found in late O dwarfs 
dMartins et al.ll2004 [2005bt iMarcolino et al.ll2009l) . 

We have fitted only the A'-band part of the observed spectra, 
but our best fit models yield the complete spectral energy dis- 
tribution from which one can compute the number of ionizing 
photons and the ionizing luminosity. The values we derived are 
listed in columns 10 and 11 of Table [4.2.21 Such quantities are 
important to understand the physics and dynamics of H 11 regions 
(see Sect. 16. U . 

4.2.2. Age determination 

Fig- 13 (left) shows the HR diagram of t he ionizing cluster 
of RC W 79. Evolutionary tracks are from iMevnet & Maeden 
(2003) and include the effects of rotation. All stars lie rather 
close to the main sequence, on a relatively narrow band, indi- 
cating that they most likely formed in a single star formation 



Cygni profiles or from the width of strong emission lines. In our 
case, none of these diagnostics are available. Consequently, we 
decided to adopt Voo. We chose a value of 2000 km s _1 as repre- 
sentative of early and mi d O stars, and 1 000 km s _I as typical of 
late O and B stars (e.g. iPrinia et al.|[l990l) . We also adopted the 
so-called /3 parameter. Our model atmospheres require an input 
velocity structure which is constructed from a quasi-static photo- 
spheric structure to which a /3 velocity law is connected. This law 
is of the form: v = v M (l —R/rf where R is the stellar radius. The 
value of 0.8 we adop ted for /3 is typical of O dwarfs/giants (e.g. 
Repo lust et al.l 2004). For the input ph otospheric structure, w e 
used the OSTAR2002 TLUSTY models dLanz & Hubeny|2003l). 



Finally, we used the solar abundances of Grevesse et al.l d2007l) 



for the elements included in our models, namely H, He, C, N, O, 
Ne, Mg, Si, S, Fe and Ni. 

Fig. IA.3I and IA.4I shows the best-fit we obtained for the 
brightest ionizing sources of the three H 11 regions. In general, 
those fits are of good quality. The only important discrepancy 
in a few stars is the Hei 2.056 yum line which is predicted too 
strong. However, this line has been shown to be ext remely sen- 
sitive to line-blanketing (N aiarro et al.l 1 1 994l 120061) . Since our 
models include only a limited number of eleme nts and since 
some atomic data for metals remain uncertain (see Naiarro et al. 
2006), the observed discrepancy is not surprising. The derived 
parameters for each star are gathered in Table 14.2^21 

Fig. [6] shows an example of mass-loss rate determination 
based on Bry in the case of star RCW79 5a. Reducing M 
produces a deeper absorption line. However, below a value of 
3 x 1 0~ 8 M yr~ 1 the line becomes little sensitive to any change of 
mass-loss rate. In practice, it is the wind density which is derived 
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Fig. 6. Mass loss estimate for star RCW79 5a. The black solid 
line is the observed spectrum. The colored lines refer to models 
with different mass loss rates. A conservative upper limit of 1CT 7 
M yr~' is derived. 



event. By finding the theoretical isochrone best representing the 
distribution of stars, the age of the population can be derived. 
For this, we use the following function: 

2 ^ (lo g (T f? -lo g{ r 7f) f ( l og( L')-log(Ln) 2 
Xis ° 2j crZ(log(P ff) ) cr\log(U)) 

where "iso" stands for "isochrone" and n is the number of stars. 
T'eff an< ^ L' so are the coordinates of the point of the isochrone 
closest to the values (T e ff, L) representing a star's position. 
Hence, x 2 iso is an estimate of the sum over all stars of the distance 
of a star's position to the selected isochrone. The isochrone for 
which this distance is minimized indicates the age of the popu- 
lation. In practice, %f is evaluated for several isochrones corre- 
sponding to ages from to 5 Myr. The results are displayed in 
Fig. [8] The minimum is found for an age of 2.0-2.5 Myr and a 
typical uncertainty of ~0.5 Myr. Note that this estimate takes all 
stars into account. If we use only the five brightest stars, we find 
a very similar result. This is because the lower part of the HRD 
is almost insensitive to age due to larger error bars and small 
separation between isochrones. 

The position of the brightest source of the compact H n re- 
gion (CHII) on the border of RCW 79 (see Figs. 4 of ZA06 and 
Fig. |B.4t is shown by the green triangle in Fig. [7] (left). It is 
essentially indistinguishable from the position of the ionizing 
sources. An age estimate for that star would give the same result 
as for the ionizing sources, within the error bars. The other stars 
of the CHII region do not help to refine this estimate since they 
are much fainter and cooler and fall in the region of the HR di- 
agram where the isochrones are tightly packed. Given the errors 
on the effective temperature and luminosity of star number 6, 
one can exclude an age younger than 0.5 Myr. This means that if 
there is any age difference between the CHII region and the ion- 



Fig. 8. Age estimate of the cluster of the ionizing sources of 
RCW 79. A value of 2.3+0.5 Myr is obtained. The solid line 
is a spline fit curve. 

izing stars of RCW 79 (as would be expected in case of triggered 
star formation) it is not larger than 2 Myr (see also ZA06). 

The right panel of Fig. [7] shows the position of the ionizing 
stars of RCW 82 and RCW 120. Only the two O stars (e2 and e3) 
are shown for RCW 82 since we have not analyzed the properties 
of the Be star el (our atmosphere code does not allow a treatment 
of circumstellar material). Unlike RCW 79, the small number of 
stars prevents an accurate age determination. Besides, the large 
uncertainty on T e fr and the relatively low luminosity of the ob- 
jects (compared to the brightest stars in RCW 79) complicates 
any attempt to constrain the stars' age. The only safe conclusion 
one can draw is that the ionizing stars are younger than 5 Myr. 
Other than that, any age smaller than 5 Myr is possible. 

From interpolation between evolutionary tracks in Fig. [7] 
one can estimate the masses of the stars of RCW 79. They are 
summarized in Table 14.2.21 Five stars have masses larger than 
30M G . If we assume a Salpeter IMF between 0.8 and 100 M Q 
for the entire population formed with the ionizing sources, one 
finds a total mass of about 2OOOM . Hence, if the star forma- 
tion event which produced the ionizing stars of RCW 79 also 
gave birth to lower mass stars, a rather massive cluster was born 
2.0-2.5 Myr ago. In that case, the faint sources seen in Fig. [4] 
are probably intermediate mass main sequence stars born at that 
time. 



5. YSOs on the borders of H n regions 

In this section we now focus on the second groups of objects: 
embedded sources on the borders of the H n regions. They 
have been classified as candidate YSOs by ZA06, ZA07, DE09 
and PO09. Our aims are the following: 1) confirm spectroscop- 
ically the YSO nature of these sources, 2) get more insight 
into their physical properties (presence of disks, envelopes, jets, 
outflows). For this, twenty-three YSOs have been observed: 8 
around RCW 79, 5 around RCW 82, and 10 around RCW 120. 
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Fig. 7. HR diagram of th e ionizing sources of RCW 79 (left) and RCW 82 and RCW 120 (right). Evolutionary tracks are from 
Meynet & Maeder (2003). Isochrones corresponding to ages of 1, 2, 3, 4, 5 Myr are shown by magenta dashed lines. Left: Blue 
circles are the ionizing stars of RCW 79 and the green triangle is star number 6 in the compact H n region on the border of RCW 79. 
Right: the ionizing source of RCW 120 is shown by the green triangle, while the two O stars ionizing RCW 82 (sources e2 and e3) 
are shown by the blue circles. The Be star (el) is not shown because its properties have not been derived in our study. 



Table 2. Derived stellar properties of the ionizing stars of RCW 79, RCW 82 and RCW 120. 



Source 



ST 



T e s 
[kK] 



[M ] 



log(Qo) 



log(Lo) 



RCW79 2 
RCW79 4a 
RCW79 5a 
RCW79 6a 
RCW79 6b 
RCW79 7 
RCW79 8a 
RCW79 8b 
RCW79 9a 
RCW79 9b 
RCW79 10 
RCW79 12b 

CHII #6 

RCW82 e2 
RCW82 e3 



07.5-8V/III 

04-6V/III 

04-6V/III 

O6-8V/III 

09-B2V 

O6-8V/III 

09-B2V/III 

09-B2V/III 

09-B2V/III 

09-B2V/III 

06.5-8V/III 

09-B2V/III 



36+2 
40±2 
40±2 
38±2 
28±4 
37±2 
30±3 
30±3 
30±3 
30±3 
38±2 
28±4 



4.95±0.22 
5.56±0.22 
5.56±0.22 
5.33±0.22 
4.15±0.27 
5.30±0.22 
4.76±0.24 
4.67±0.24 
4.57+0.24 
4.53±0.24 
5.26±0.21 
4.21±0.27 



26.2^T 
46 1 +133 

-t?j 

36.2$ 
13.2!| 
35.1+S| 
19.2« : ° 
18.6!« 
17. 7 +43 
17.3^^ 
34 5 + <^ 
36^ 



48.45+0.23 
49.25±0.23 
49.25±0.23 
48.91 ±0.23 
47.00±0.31 
48.80±0.23 
47.33±0.26 
46.75±0.26 
47.95±0.26 
48.52±0.26 
48.96±0.22 
47.06+0.31 



09-B2V/III 
09-B2V/III 



6.8 

34±2 4.95±0.27 24A + _f Q 
35±2 5.07±0.27 27.5 + 



RCW120e O6-8V/III 



37.5±2 5.07±0.21 



-5.8 



29.9 + " 



37.90±0.23 
38.74±0.23 
38.74±0.23 
38.40±0.23 
36.42±0.31 
38.28±0.23 
36.75±0.26 
36.16±0.26 
37.38±0.26 
37.99±0.26 
38.45±0.22 
36.48±0.31 



07.5-9.5V/III 37.5±2 5.25±0.22 34.1+" 48.79+0.23 38.26±0.23 



48.37±0.28 37.81±0.28 
48.59+0.28 38.05+0.28 

48.58+0.22 38.05+0.22 



Note 1: The columns are: target name, spectral type, effective temperature, luminosity, evolutionary mass, ionizing flux, ionizing luminosity. 
Note 2: The errors on the last two parameters only take into account the uncertainties on T c ff and log j^. 



They were selected mainly from their high A'-band magnitude. 
Their position and photometry is given in Table [3] Appendix 151 
gives further details on each individual objects. We first present 
a general description of their spectral properties (Sect. I5.lt be- 
fore confronting them to near and mid infrared photometry (Sect. 
15. 2t . We then focus on YSOs with strong H2 emission lines in 
Sect. 15.31 The morphological and kinematic properties of a few 
objects are presented in Sect. 15. 41 where we show that the YSOs 
are associated to the HII regions. 



5. 1 . Spectroscopy of embedded sources 



The 23 objects we have observed on the borders of the three H 11 
regions show four main type of spectral signatures: Bry emis- 
sion, H2 emission, CO bandheads emission and CO bandheads 
absorption. Some sources are entirely dominated by one type of 
spectroscopic signatures, but others can show several features, 
such as Bry and H2 emission, or H2, Bry and CO emission. Table 
[4] summarizes the main spectroscopic properties of each source. 
Their spectra are shown in Figs. IA.5I and IA.6I Below, we de- 
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scribe the various sources, grouping them according to the main 
feature they display. 

• H2 emission sources: ro-vibrational H2 transitions are ob- 
served in these objects. The strongest one is H2 1-0S(1) 
at 2.122 yum and it is sometimes the only line observed. 
However, in the majority of "H2 emission sources" a number 
of other lines are observed. They are identified on Fig. IA.5I 
(top right). Such lin es are known to be formed either ther- 
mally in shocks (e.g. lBurton et al.|[l990h . or by fluorescence 
after excitation of H2 molecules by n on ionizing FUV radia- 
tion i n the Lyman- Werner band (e.g. lBlack & van Dishoeckl 
119871) . In the latter case, the lines are formed under non LTE 
conditions. They can thus trace either outflows or photo- 
dissociation regions (PDRs) in young star forming regions. 
In Sect. 15.31 we will analyze quantitatively these sources to 
track the nature of the exciting mechanism. 

• Bry emission sources: strong Bry emission is observed in 
several sources, sometimes in combination with H2 emis- 
sion. Bry may have various origins. According to Bik et 
al. (2006), the large line widths (100-200 km s" 1 ) could 
be explained by ionized gas flowing from the surface of a 
circumstellar rotating disk, rather than a pure nebular ori- 
gin. Other mechanisms such as accretion flows or outflows 
are alternative possibilit ies. In the form er case, magneto- 
spheric accretion ([v an den Anckerj 120051) or inner disk ac- 
cretion dMuzerolle et al. 20041) theoretical models predict 
Bry emission. Fo r outflows, classi cal stellar winds, X-winds 
dShuet alj [19941) or disk winds dBlandford & Pavnd [1981 
e.g.) lead to ionized hydrogen emission. Interferome t ric ob- 
servations of Herbig Ae/Be stars by iRraus et al I (120081) 
favour Bry emission in extended stellar or disk winds. 

• CO bandheads emission sources: CO first overtone emis- 
sion at 2.2935 fxm (CO(2-l)), 2.3227 //m (CO(3-l)), 2.3535 
fim (CO(4-2)), 2.3829 fim (CO(5-3)) is observed in eight 
sources. This spectral signature has been observed in a 
number of low and high mass YSOs dScoville et all [1971 
IChandleretaTJ 119931 iBik et alj 12006b . It is commonly at- 
tributed to_Jrradjationof_an accretion disk by the central 
object dChandler et al.lll995l: Uraus et alj|2000b . This inter- 
pretation is motivated by the temperature (2500 to 5000 K) 



and density (n ~ 10 cm" 3 ) conditions required to pro- 
duce such an emission, placing the emitting zone close to 
the ce ntral star (within 0.5 AU for the Be star 5 1 Oph stud- 
ied bv lBerthoud et alj|2007l [Tatulli et alj|2008l) . The geom- 
etry of the bandhe ads depends on the Keplerian rotation 
and viewing angle . IBik & Thl d2004l) used the formalism of 
IRraus et aiTd2000l) to show that CO emission in the massive 
YSOs observed bv lBik etaU d2005l) was produced in a kep- 
lerian disk. Our spectra do not have high enough resolution 
to detect line asymmetries. In the strongest CO emitters of 
our sample, Nai 2.206-2.209 fjm and Cai 2.261-2.264 yum 
emission is also detected (see Fig. IA.51 ). 
CO bandheads absorption sources: the CO first overtone 
lines are observed in absorption in a few sources (Fig. |A.5l 
bottom right). The line depth ranges from less than 10% to 
about 40% of the continuum.. Nai 2.206-2.209 //m, Cai 
2.261-2.264 /mi and Mgi 2.281 fim lines are usually ob- 
served in absorption as well. In one object (RCW120 51b) 
a significant Bry emission is observed, while in another one 
(RCW82 98) H2 emission lines are present. CO absorption in 
the Zf-band is typical of cool stars (red supergiants, AGBs, 
giants), but it is also detected in YSOs (jC asali & Ma tthews! 
ll992tlHofrrneister et alj2 006: Aspin et alJ2009al) . In the lat- 



ter case, it is thought that CO absorption is observed due to 
the absence of a strong disk and/or envelope emission which 
otherwise dominates the K-band spectrum and leads to ei- 
ther featureless or emission line spectra. This is consistent 
with the fact that CO ab sorption YSOs tend to b e class II 
(and not class I) objects dCasali & Matthews|[l992l) . CO ab- 
sorption is also observ ed in the rare class of e ruptive FU 
Orionis objects (FUOr. iReipurth & Aspin|[l997l) . They are 
thought to be active T-Tauri stars experiencing an accre- 
tion burst. lAspin et all d2009bl) collected spectra of VI 647 
Orionis, the FU Or star responsible for the appearance of the 
McNeil nebula dMcNeil et alj|2004l) . and show that its CO 
bandheads shift from abs orption to emissio n when the star 
experiences a FUOr event. lCalvet et al.1 dl99ll) explained the- 
oretically this type of transition by an increase of the accre- 
tion rate. EXor variables which are pre main sequence ob- 
jects similar to Fuor stars but wi th weaker outbursts als o dis- 
play CO absorption bandheads dLorenzet ti et al. 2009). 

All these spectral signatures are typical of YSOs. In spite of 
this general grouping in four categories, we stress again that a 
few sources share several of these spectral characteristics, high- 
lighting the complexity of the physical phenomena associated 
with YSOs. In the following, we attempt to better characterize 
the nature of the sources using infrared photometry. 



5.2. Near and mid infrared photometry 

All our sources have been observed by 2MASS. We have seen 
that in the case of the ionizing stars of RCW 79, some 2MASS 
sources were actually multiple. Similarly, some YSO candi- 
dates have been resolved into several sub-components by our 
SINFONI observations (see Figs. [OE2] and ED. Hence, the 
2MASS photometry is only indicative for those specific sources. 
Fig.[9]shows a JHK color-color diagram based on 2MASS pho- 
tometry. Star symbols refers to multiple sources. When avail- 
able, we have used ESO/NTT/Sofl photometry from ZA06, 
ZA07 and PO09 since it is more accurate and does not suffer 
from the same uncertainties. The main sequence is also plotted 
together with Ay=40 mag reddening vectors. The main conclu- 
sion is that most of the objects show near-IR excess. Eight ob- 
jects lie close to the reddening vectors. Out of them, three are 
multiple sources. Among the five remaining objects, RCW120 
C4-64 shows CO absorption bandheads and is most likely a fore- 
ground giant star. The near-IR excess observed in the objects far 
from the reddened main sequence is usually indicative of the 
presence of an envelope or a disk. It thus confirms that most of 
our sources are YSOs. 

To get more insight into the nature of the sources, we have 
retrieved Spitzer/IR AC photometry (when available) from the 
GLIMPSE survey dBeniamin et all 120031) to build the color- 
color diagram shown in Fig. [10] The different symbols corre- 
spond to the four spectroscopic groups defined in Sect. 15. II Class 
I and II objects gather within the dashed and solid lines respec- 
tively. Stars are shown by pentagons close to (0,0). 

All our sources fall in the area corresponding to Class I 
and Class II YSOs. They gather around the color point (0.8, 
0.8); some of them may be Class II so urces displaced by ex- 
tinction (see Fig. [TOj. As discussed by iRobitaille et alj ([2008 ) 
a few of these sourc es may be extreme AGB stars (see also 
ISrinivasan et al.ll2~009l) . It is probably the case of the CO absorp- 
tion sources RCW 120 39 and 64. Four sources have different 
colors: two are pure H2 sources (RCW79 237 and 3C3), the two 
others have colors typical of PDRs (RCW79 1C4, associated 
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Fig. 9. JHK Color-color diagram (magnitudes from 2MASS and 
NTT-SofI) of the YSOs associated with RCW 79, RCW 82 and 
RCW 120. The location of main sequence stars is shown by 
the solid line. Reddening lines have a length corresponding to 
a visual extinction of 40 mag. Filled (open) triangles correspond 
to good (uncertain, from 2MASS catalog flag) magnitude mea- 
surements. Star symbols following the names indicate multiple 
sources. The typical error bar of good measurements is indi- 
cated. 



with extended 8 /urn and 24 fim emission, and RCW79 2C3). 
The two sources with the largest [3.6]-[4.5] co lor are the ones 
with th e strongest H2 lines (Fig. [A3] top right). Smit h & Rosenl 
(2005) showed that numerous H2 lines contributed to the IRAC 
fluxes, but that the integrated line intensity was the largest in 
the 4.5 fim band. Hence, the presence of strong H2 lines in the 
K-bmd of RCW79 3C3 and RCW79 237 is probably respon- 
sible for an intense emission in the 4.5 /im band, leading to an 
increased flux compared to the 3.6 /im channel. Consequently, 
the value of [3.6]-[4.5] is larger in those two objects, explain- 
ing their location in Fig. [TO] The other H2 dominated sources 
have weaker lines, and consequently their IRAC colors are not 
affected. 

Except for the four cases discussed above, the general con- 
clusion is that objects with similar IRAC colors can present dif- 
ferent spectroscopic signatures. Said differently, there is no di- 
rect correlation between IRAC colors and a spectroscopic group. 
Once again, the multiplicity of a few sources complicates the in- 
terpretation. 

Fig. [TT] shows a color-ma gnitude diagram based on 
Spitzer/MIPS 24 fim photometry dCarev et al] |2009). comple- 
mented by AKARI observations (Zavagno et al. in prep). The 
magnitude in this band corrected for the distance modulus of 
each source (noted [24]^), is plotted as a function of K— [241 . 
For comparison, the theoretical SEDs of iMartins etaT] ((2005a) 
have been used to calculate the positions of O stars in this di- 
agram. They are shown by the pentagon symbols. A reddening 
ve ctor corresp onding to Ay=40 mag and to the extinction law 
of lLutzl d 19991) is indicated. All YSOs are brighter at 24 jjm and 



Fig. 10. Color-color diagram with SPITZER/IRAC magnitudes 
of the YSOs associated with RCW 79, RCW 82 and RCW 120. 
Different symbols are used for the four spectroscopic groups de- 
fined in Sect. 15.11 triangles for H2 sources, squares for CO ab- 
sorption sources, circles for CO emission sources and diamonds 
for Bry dominated sources. Pentagons indicate main sequence 
massive stars (from Martins et al. 2005a). Colors of class II (re- 
spectively I) sources are deline ated by solid (resp. dashed) lines 
according to lAllen et al.l (120041) . The reddening vector shown by 
the black arrow is comp uted for a visual e xtinct ion of 40 mag- 
nitudes with the law of Ind ebetouw et alj J2005). Star symbols 
following source names indicate multiple sources. 



much redder (in A"-[24]) than O stars. This indicates once more 
that they have a strong infrared excess due to envelopes and/or 
disks which dominate the SED at those wavelengths. The H2 
sources seem to stand out in this diagram. While most sources 
are grouped in a region defined by -12 < [24]d c < -6 and 7 < 
^T-[24] < 10, four out of the five H2 sources for which infrared 
photometry is available are much brighter and redder ([24]rf c < 
-12 and ZT-[24] > 10). This might be an indication that those 
objects have larger amount of circumstellar material tracing an 
earlier state of evolution, or that they have larger masses. 

This is confirmed by Fig. Q~2] showing the [3.6]-[5.8] vs 
[8.0] -[24] diagram. Here again, the H2 objects are clearly lo- 
cated towards redder colors (with the exception of RCW120 76). 
The solid lines delinea tes the location of stage I and II objects of 
Robitaille et al.l d2006l) . This classification is the theoretical ana- 
log of class I and II sources, meaning that stage I objects corre- 
spond to models with an envelope, while stage II sources have 
disks and possibly a tenuous envelope. H2 sources are therefore 
consistent with being stage I objects as suspected previously. 
The bulk of the other objects are again located in the transition 
region between stage I and stage II sources. Only RCW 120 51 
and RCW82 51 might appear as stage I objects. We will dis- 
cuss further these results in Sect. 16.21 using additional findings 
presented below. 

The main conclusion of this analysis of color-color and 
color-magnitude diagrams is that the objects dominated by H2 
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Fig. 11. Color-magnitude diagram with SPITZER/MIPS 24 fj.m 
magnitudes of the YSOs associated with RCW 79, RCW 82 and 
RCW 120. On the Y axis, the magnitudes have been corrected 
for distance. The symbols have the same meaning as in Fig.lTOl 
The arrow indicates an extinction of Ay =40 magnitudes (law of 
lLutj!999[). The pentago ns are O stars from the theoretical SEDs 
of iMartins et al.l (l2005al) . Star symbols following source names 
indicate multiple sources. 



emission seem to be different from the remaining of the sources 
we observed. They probably have an optically thick envelope. 
The other objects are more consistent with YSO surrounded by 
disks. Another conclusion regarding these sources is that they 
can have different spectroscopic signatures but similar colors. 
This can be partly explained by geometrical effects. Indeed, de- 
pending on the inclination and size of the putative disks, different 
regions are visible, and thus different spectral line s are expected. 
The corre sponding photometry is also modified. Whitne v"et al.l 
d2003blah showed that in the [3.6]-[4.5] vs [5.8]-[8.0], a spread 
of about one magnitude can be expected due to inclination ef- 
fects. This corresponds to what we observe. Hence, among the 
CO emitters, absorbers and Bry emitters, one probably sees a 
population of class II objects at different inclinations. 

5.3. Nature of H 2 sources 

Several sources observed on the border of the H n regions show 
H2-dominated spectra (see top right panel of Fig. [A3] for line 
identification). As described in Sect. 15.11 H2 emission is due ei- 
ther to thermal processes or to fluorescence. The relative inten- 
sity of individual lines is different in either cases. Comparing H2 
line strength can thus lead to the underlying excitation mech- 
anism. For a thorough comparison, line fluxes must be dered- 
dened, which requires the knowledge of extinction. This piece 
of information is difficult to obtain since the intrinsic SED of 
our sources is not known. From Fig.[9]one can see that most ob- 
jects have a visual extinction in the range 10-40 mag. We have 
thus adopted a typical extinction A^ = 2.5+1.5 for the present 



CO 
LD 

CO 

co 



- 




[8.0]-[24] 

Fig. 12. [3.6]-[5.8] vs [8.0]-[24] diagram. The symbols have 
the same meaning as in Fig. \W\ The arrow indicates an extinc- 
tion of Ay=40 magnitudes. T he solid lines de l ineate s the areas 
where stage I and II objects of Robitai fTe et al.l (120061) are found 
(see their Fig. 23). Star symbols following source names indicate 
multiple sources. 



study. We used the Galactic extinction law oc A 1 1 similar to 
Mon eti et al.l d200lh to deredden the line fluxes. 

We have computed the column densities N; of several ex- 
cited levels and plotted them as a function of level energy 
This has been shown to be a powerf ul tool to distinguish be tween 
thermal and radiative emission (e.g. lHasegawa et aljl987l) . If H 2 
emission is thermal, the values of N\jg; {gj being the statistical 
weight) follow a Boltzmann distribution and are thus aligned in a 
log(Nj/gi) - Ej diagram. Alternatively, radiative pumping should 
increase the population of high energy levels and introduce a 
departure from linearity in the aforementioned diagram. We 
have computed the column densities o f several levels using the 
same formalism and molecular data as Ma rtin-Hernandez et al.l 
d2008l) . Our results are shown in Figs. [13] The solid and dotted 
lines correspond to a linear fit of the v = 1 levels. The slope 
is directly related to the temperature of the Boltzmann distribu- 
tion best representing these energy levels. Solid lines are for the 
case Ak-2.5, while the dotted lines are for Ak=1.0 and A^=4.0. 
When extinction varies, line fluxes are affected, which explains 
the different slopes for different A K . From Fig.[T3lRCW79 3C3, 
RCW120 103 and RCW79 91b appear consistent with pure ther- 
mal emission. For RCW 79 2C3, RCW120 76, RCW120 49 and 
RCW120 107, a marginal deviation from pure thermal emission 
is found, although the bulk of the emission seems to be thermal. 
Source RCW82 7 shows too few lines to draw any final conclu- 
sion. Besides, we will see in Sect. |5.4| that H2 emission is much 
extended and patchy. Finally, RCW79 237 clearly shows a sig- 
nificant non thermal contribution. 

Two more sources (RCW79 89 and RCW79 91) display H 2 
lines, but not enough to build excitation diagrams. For these 

sources, we calculated the ratios j^jj- and 0i (defined by 
Smith et a i] ll997h . Under LTE conditions, the first ratio should 
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Table 3. 2MASS and NTT/SofI JHK, Spitzer-GLIMPSE, Spitzer-MIPS and AKARI magnitudes of YSOs. 



Star 


RA 


DEC 


J H K s [3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] 


nature 




h m s 
















YSOs around RCW 79 



237 


13:41:02.3 


-61:44:16 


17.41 


15.79 


14.23 


11.02 


9.00 


7.52 


6.49 


1.61" 


Class I 


91 


13:40:57.6 


-61:45:44 


16.62 


14.45 


12.59 


8.66* 


7.60* 


7.16 


6.21 


-0.84" 


Class I 


18 


13:40:58.3 


-61:46:50 


17.11 


13.83 


11.45 


8.57 


7.54 


6.70 


5.87 


2.21 


Class I 


89 


13:40:59.0 


-61:46:51 


18.58 


15.35 


13.34 


10.83 


10.71 










186 


13:40:58.5 


-61:46:48 


18.96 


15.86 


13.91 














2C3 


13:40:32.7 


-61:47:20 


14.08 


12.79 


11.82 


9.54 


9.04 


6.62 


4.96 


-0.87" 


Class I 


3C3 


13:40:26.6 


-61:47:56 


17.24 


14.32 


12.57 


9.21 


7.02 


5.84 


5.30 


0.26" 


Class I 


1C4 


13:39:53.6 


-61:41:21 


12.19 


11.45 


11.01 


8.19 6 


7.94* 


5.09* 


3.41* 


0.26" 




4C4 


13:39:55.9 


-61:40:58 


14.68 


13.54 


10.62 


7.27 


6.22* 


5.28 


4.67 


2.01* 


Class I 


6C4 


13:40:08.9 


-61:41:38 


13.91 


11.96 


10.29 


8.78 


8.18 


7.51 


6.38 


2.21" 


Class I 


YSOs around RCW 82 


4 


13:59:34.8 


-61:19:22 


17.59: 


15.06 


12.41 


9.89 


9.15 


8.59 


8.30 


5.04 


Class I 


7 


13:59:47.3 


-61:22:02 


18.77 c 


15.85 c 


14.03 c 


11.34 


10.14 


9.02 


8.17 


4.05 


Class I 


35 


13:59:50.9 


-61:25:22 


14.29 


12.51 


10.96 


8.61 


7.79 


6.97 


6.11 


2.58 


Class I 


51 


13:59:57.6 


-61:24:37 


17.81 c 


13.80 c 


11.13 r 


8.31 


7.48 


6.01 


4.80 


-0.19 


Class I 


98 


13:59:36.7 


-61:20:48 


16.27: 


15.31: 


13.57 


10.09 


9.24 


8.47 




4.01 




YSOs around RCW 120 


39 


17:12:00.9 


-38:31:40 


17.29: 


14.95: 


11.91 


7.43 


6.57 


5.66 


4.93 


2.26 


Class II 


49 


17:12:08.0 


-38:30:52 


16.24 


13.67 


12.18 


10.59 


9.53 


8.85 


8.16 


4.36 


Class I/II 


51 


17:12:28.5 


-38:30:49 


12.85 


11.73 


10.97 


10.00 


9.20 


8.31 


7.78 


2.88 


Class I/II 


C4-64 


17:12:40.8 


-38:27:29 


13.37 


11.85 


11.06 


9.78 


9.17 


8.50 


7.60 


3.75 


Class I/II 


C4-67 


17:12:41.1 


-38:27:17 


13.03 


10.96 


9.51 


7.78 


7.03 


6.27 


5.56 


2.09 


Class I/II 


C4-69 


17:12:42.2 


-38:26:58 


15.16: 


13.50 


12.28 


10.82 


9.75 


8.78 




3.87 


Class I/II 


76 


17:12:46.0 


-38:25:25 


17.03: 


13.66 


10.82 


7.61 


6.40 


5.37 


4.47 


1.37 


Class I/II 


103 


17:12:40.2 


-38:20:31 


14.23 


12.57 


11.31 


9.40 


8.29 


7.24 


6.13 


2.87 


Class I/II 


106 


17:11:50.7 


-38:19:55 


17.62: 


15.80: 


12.41 


8.27 


6.89 


5.79 


4.95 


0.97 


Class I 


107 


17:12:33.6 


-38:19:52 


16.03: 


12.85 


10.19 


7.13 




5.32 


4.57 


1.48 


Class I/II 



Note: Classification into class I/II objects is from ZA06, PO09, DE09. The notation ":" uncertain flux measurements due to the weakness of the source. 
The typical errors are 0.1 mag (respectively 0.3 mag) for IRAC (resp. MIPS) photometry. 
" AKARI magnitude 

* Magnitude obtained by aperture photometry (see ZAV05) 
NTT magnitude 



have values between 5 and 100 (for te mperatures ranging be - 
tween 3000 and 1000 K respectively, see lGustafsson et al.l2 008). 
and the second ratio should be close to 3. When fluorescence 
dominates, lower values are expected. For RCW79 89, the ob- 
served values are 4.61+0.79 and 2.89+0.56 respectively, indi- 
cating that H2 has probably multiple origins. For RCW79 91, 
FUV excitation is clearly favored since !' = 6.33+0.14 and 

*2-lS(l) 

(pi = 0.99+0.14. 

In conclusion, most sources show a thermal H2 spectrum 
with, in some cases, a possible contribution from fluorescence. 
Three sources (RCW79 237, RCW79 89 and RCW79 91) are 
probably dominated by non thermal emission, most likely due 
to the presence of n earby PDRs. Therm al emission is usually at- 
tributed to shocks. iDavies et a l. (2000) summarized the various 
types of sho c ks enc ountered in the ISM using the models of 
iBurton et alj (|T990). They note in particular that fast J-shocks 
are usually dissociative. In that case, molecules are formed on 
grains in the post-shock region and are produced in excited 
states. Consequently, their emission spectrum is more typical 
of a non-thermal mechanism since the molecules will cascade 
down to lower energy levels. Hence, fast J-shocks seem to be 
excluded to explain the thermal emission of most of our H2 
sources. Similarly, slow C-shocks are not favored either since 
they are associated with rather low temperatures (<300 K). From 
Fig [13] we see that temperatures of 1000-2500 K are estimated, 

3 Two types of shocks are usually encountered: J-shocks (for Jump 
shocks) in which the change in density and velocity occurs on scales 
much shorter than the radiative cooling scale, and C-shocks (for 
Continuous shocks) in which the changes take place on a longer scale. 



inconsistent with slow C-shocks. Hence, slow J-shocks or fast 
C-shocks, both heating the gas to 2000-3000 K, are preferred to 
account for the observed thermal H2 emission in our sources. 



5.4. Morphology and kinematics of selected YSOs 

In this section we gather information on the morphology and 
kinematics of some YSOs in order to study their physical as- 
sociation with the H 11 regions and to get more insight into the 
origin of the observed emission lines. 

We first focus on the kinematics of the YSOs themselves. 
The emission lines of their spectra (Bry and/or H2) have been 
fitted by ID Gaussians using QFitsView. We have first measured 
the full width at half maximum (FWHM) of the sky lines to 
estimate the spectral resolution of the observations. We obtain 
FWHM= 84 + 6 km s -1 , thus a resolution R= 3575 ± 250. Then 
we have measured the central positions and the widths of the H2 
and Bry lines (only for the lines presenting a good signal to noise 
ratio). These quantities are given in Table[5] The first and second 
columns identify the YSO; the indication "ext" means that we 
have measured the nebulous extended emission in the vicinity of 
the YSO. The following columns give, for a selection of lines, 
the LSR velocity (Vlsr) and the intrinsic line width (AV). The 
absence of measurement for some lines is due to the poor qual- 
ity of the data at hand. We estimate that the uncertainty on Vlsr 
is ~ 10 km s _1 . The relative uncertainty on the measured line 
widths if of 10%, but the uncertainty on AV is large. Its is about 
20 (40) km s _1 for wide (narrow) lines. This is mainly due to the 



Table 4. Characteristics of the YSO spectra. When present, lines are in emission unless stated differently. 



Source 


Bry 


H 2 


CO 


Notes 


YSOs around RCW 79 


237 


no 


2.0338, 2.1218, 2.2233, 2.2477 


no 




237-b 


no 


no 


absorption 


Na i, Ca i absorption 


91 


yes 


2.0338, 2.1218, 2.2233 


no 




91-a 


yes 


2.033,2.1509 


emission 




91-b 


yes 


1.9576, 2.0338, 2.1218, 2.2233 


weak emission 




18 


yes 


no 


emission 


Na i emission 


186 


yes 


no 


no 


+3198 


89 


yes, weak 


no 


no 




2C3 


yes, on top of absorption 


2.0338, 2.1218, 2.2233 


no 




3C3 


yes, weak 


2.0338, 2.1218, 2.2233, 2.2477 


no 




3C3-b,c,d 


no 


2.1218 


no 




1C4 


yes, absorption 


no 


absorption 


Nai, Cai absorption. PDR associated with this objet 


4C4 


yes, weak 


no 


no 


faint companions 


6C4 


yes 


no 


no 




YSOs around RCW 82 


4 


yes 


no 


no 




7 


? 


1.9576, 2.0338, 2.1218, 2.2233 


emission 




35 


yes 


2.0338, 2.1218, 2.2233 


emission 


Na i emission 


51 


yes 


2.0338, 2.1218 


emission (weak) 




98 


no 


1.9576, 2.0338, 2.1218, 2.2233, 2.4066, 2.4134 


absorption 


H2 lines probably not related to YSO, see Fig.IB.ll 


YSOs around RCW 120 


39 


no 


no 


absorption 


Na 1, Ca 1 absorption 


49 


yes 


2.0338, 2.1218, 2.2233, 2.2477 


emission 




49-a 


no 


2.0338, 2.1218, 2.2233 


absorption 




49-b 


no 


2.0338, 2.1218, 2.2233 


absorption? 




51-a 


yes 


no 


no 




51-b 


no 


no 


absorption 




51-c 


no 


no 


absorption 


Na 1, Ca 1 absorption 


C4-64 


absorption 


no 


absorption 


Na 1, Ca 1 absorption 


C4-67 


yes 


2.1218 


emission 




C4-69a 


no 


no 


no 


featureless 


C4-69b 


yes 


no 


absorption 




C4-69c 


no 


no 


absorption 


Na 1, Ca 1 absorption 


C4 A 


absorption 


no 


absorption 


PDR associated with this object 


C4Ba 


absorption 


no 


no 


PDR associated with this objet 


C4Bb 


no 


no 


absorption 


Na 1, Ca 1 absorption 


76 


no 


2.0338, 2.1218, 2.2233 


no 




103 


yes 


1.9576, 2.0338, 2.1218, 2.2233, 2.4066, 2.4237 


emission 


Nai, Cai emission 


106 


yes 


no 


no 




107 


yes, weak 


2.0338, 2.1218, 2.2233 


no 





14 



F. Martins et al.: Near-IR integral field spectroscopy of ionizing stars and young stellar objects on the borders of H n regions 



i i i i | i i i i | i i i i | i i i i- 
jfn' RCW79 3C3 I 
fir? 



^.3 . 




2025 



\\f\- -2 - 
k •• V : : l 



I I "l I I V I I I H 





■1 



-1 1 \S1 1 M II 1 


i i i i i i i i i i- 


I \ RCW79 237 I 






I 1062 k'v 




-i 1 1 1 1 1 1 1 1 i'A Ki 1 1 1 1 1 1 1 1- 



- 

-1 r 
-2 - 



I — I I -JV.I I I I I I I I I I I I I I I I I I I- 
RCW79 2C3 : 




1197 K •, 

H l I l l l l I l l i- l\l l 'l l l I l l l l- 



0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 



GO 

\ 

QjO 
O 



l l l l I l l l l | l l l l | l l l l- 
RCW120 76 I 




I I I I I I I I I I I 11 I M I I I I H 



I.J I I I I | I I I I | I I I I | I I I I- 

RCW 120 103: 




2765 K 
1 1 1 1 1 1 1 1 1 1 1 1 i 1 1 f.l i 1 1 i- 



M'L^.I I I I I I I I I I I I I I I I I I I. 

RCW120 49 : 




-2 - 

I 1573 K •• 



I I I I I I I I I I I M M I I -I I I I- 



0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 



rjvj. g i i | i i i i | 1 1 i i | i i i i-i 
RCW82 7 




- -1 - 



•. - -2 - 



-1 1 111 1 


III 


II Ml 1 1 Ml II- 




RCW120 107 1 






1 

? j 


I 945 


K '\ 




^ 1 1 1 1 1 


■ III- 


1 1 1 1 1 1 1 1 1- 



|-l 1VJ. I I I I | I I I I | I I I I | I I I I- 

RCW79 91b : 




i- -i i\ 1 1- - 



0.5 1 1.5 2 2.5 



0.5 1 1.5 2 2.5 
Ej [10 4 K] 



0.5 1 1.5 2 2.5 



Fig. 13. H2 excitation diagram of the sources showing H2 emission spectra (either pure or mixed with other spectral features). 
All N/g numbers are relative to N/g for transition 1-0S(1). The solid line are fits to the v = 1 level populations assuming LTE 
and A K =2.5. The corresponding temperature is indicated for each YSO. The dotted lines are fits for the cases where A^=1.0 and 
Aa:=4.0. The error bars take into account measurement errors and an uncertainty of +1.5 in Aa:. 



large uncertainty on the instrumental width. Hence, values given 
in Tableware only indicative. 

The Vlsr velocities can be compared to the Vlsr velocities 
of the ionized gas and associated molecular material. These ve- 
locities are in the ranges -52/-46 km s" 1 for RCW 79 (ZA06), 
-52/-42 km s" 1 for RCW 82 (PO09), and -14/-8 km s" 1 for 
RCW 120 (ZA07). The comparison shows that (and taking into 
account the uncertainty on the velocity): 

• the YSOs velocities in RCW 79 are in rather good agreement 
with those of the H 11 region, except for 3C3. In this YSO 
the H2 gas approaches us with a velocity of the order of 30 
km s _1 . 3C3 has a nearby extended structure (a bow-arc like 
structure, see Fig. IB. 11 ) presenting the same velocity. 

• the YSOs in RCW 82 have velocities comparable to those 
of the associated region; an exception is source 35 with very 
different velocities for the H2 and Bry lines. The H2 material 



approaches us with a velocity ~50 km s with respect to the 
material emitting the Bry line and the RCW 82 region. 
• in RCW 120, YSO #51 has a velocity comparable to that of 
the H 11 region; it is the only one (see also below). YSOs #49, 
#76, and #107 are approaching us with a velocity of about 20 
km s . The situation is worse for YSO #103, and resembles 
very much that of source 35 in RCW 82: the Bry line has 
the same velocity as the RCW 120 region, whereas all the 
H2 lines show an approach velocity of some 50-60 km s . 
This may be indicative of flows (see below), but may also 
suggest that these YSOs are not associated with RCW 120. 
We do not favor this last solution as these three YSOs are 
observed in the direction of well defined condensations in 
the shell of dense neutral material interacting with the ion- 
ized region. Molecular observations of this shell (allowing 
velocity measurements) should help to clarify this point. 
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Table 5. Velocity measurements. For each line, the first value is the LSR radial velocity while the second value (in parenthesis) is 
the line width. 



Region 


source 


Bry 


H 2 1.9576 


H 2 2.0338 


H 2 2.1218 


H 2 2.2233 


H 2 2.4066 


H 2 2.4237 


RCW79 


18 


-58 (196) 








— 




— 




91 


-75 (202) 


-- 


-- 


-57 (60) 


-56 (-) 


— 


-- 




91a 


-63 (152) 




-- 






-- 


-- 




91b 




— 




-58 (72) 


-56 (80) 


— 


— 




237 ext 


-57 (-) 


-- 




-53 (-) 


— 




— 




237 


— 


-- 


-62 (83) 


-60 (78) 


-62 (-) 


-59 (63) 


-56 (75) 




1C4 


-49 (-) 


-- 


— 


-46 (-) 


— 


— 


— 




2C3 




-- 


-43 (-) 


-45 (33) 


-46 (nr) 


-44 (nr) 


-45 (nr) 




3C3 


-- 


-95 (95) 


-80 (103) 


-79 (88) 


-76 (75) 


-79 (84) 


-77 (76) 




3C3 bow 




-88 (-) 


-73 (-) 


-72 (72) 


-71 (70) 


-67 (50) 


-69 (49) 




CHII 


-42 (-) 






-49 (-) 








RCW 82 


4 


-52 (228) 








— 








7 ext 


— 


— 


— 


-56 (43) 


— 


— 


— 




35 


-50 (209) 


— 


-- 


-100 (188) 


— 


— 


— 




51 ext 








-56 (nr) 










51 


-74(181) 






-59 (66) 










98 






-57 (-) 


-55 (69) 


-56 (-) 


-47 (45) 


-51 (nr) 


RCW 120 


49 


-29(162) 






-30 (142) 










51 ext 


-10(-) 
















51 


-5 (162) 
















76 






-30 (-) 


-28 (78) 










103 


-12(137) 


-67 (-) 


-74 (-) 


-79(111) 


-76 (-) 








107 


-30(171) 






-(103) 









Note: All velocities are in km s The notation "nr" stands for "not resolved" and "ext" refers to extended emission. 



Table [5] shows another interesting feature: the Bry emission 
lines are wide (up to 200 km s -1 ), and always wider than the H2 
lines. This again probably indicates that they do not originate 
from the same zones of the YSOs. 

These results indicate that the geometry of the YSO is rather 
complex, at least in some cases. One can gain further insight into 
those properties by looking at emission maps and velocity maps. 
Fig. [14] to [16] show some illustrations. In those figures, the con- 
tours of the integrated A"-band emission are shown on top of im- 
ages centered on Bry and/or Fb 2. 12 yumemission lines (and from 
which the continuum is subtracted). Velocity maps built from H2 
2.12 fim are also shown when they could be extracted from the 
data. The main conclusion regarding morphology is twofold: 1) 
Bry emission is always centered on the emission peak of the in- 
tegrated A"-band image, i.e. it coincides with the position of the 
YSO; 2) H2 2.12 /mi is emitted on a much wider scale, with 
emission peaks sometimes not corresponding to the position of 
the YSO. Three YSOs have been selected to illustrate these re- 
sults, but they can be generalized to all YSOs (see appendixIBl. 

The case of RCW 120 51 is spectacular. The main A"-band 
source in this region is 51a (see Fig. [T4l. Its spectrum is domi- 
nated by Bry (Fig. IA.5t and Fig. [T4l confirms that is is the main 
emitter in this line. However, we note the presence of a Bry 
"cone" pointing towards 51a and, more generally, towards the 
center of the Ft n region. This structure is very typical of the ion- 
ization front in the vicinity of RCW 120 51. Measurements of 
the radial velocity along this cone (~ -6 km s _1 ) indicates values 
similar to that of 51a and to that of the H 11 region. Hence the 
YSO and the Bry cone are physically associated, and are also 
adjacent to the ionization front. The H2 map of Fig. [141 shows 
the existence of another emission cone, the orientation of which 
is the same as that of the Bry emission cone. However its top is 
not located on source 51a, but rather on a faint A"-band source 
located south and east of it (see Fig.[T4l. The H2 velocities along 
this cone are also consistent with those of the H n region. We are 



thus directly probing the structure of the ionization front: seen 
from the H 11 region, the line of sight crosses first a Bry cone 
tracing the ionized gas behind which is located the H2 cone lo- 
cated closer to the PDR and the molecular material. Along this 
H2 cone structure, we measure a ratio ii-2!ii (respectively 
of 6.2+0.1 (respectively 12.0+0.1) intermediate between pure 
thermal emission and radiative excitation. The presence of non- 
thermal emission is consistent with the proximity to the PDR. 
All in all, we have thus a direct proof that the YSO RCW 120 
51a is located on the ionization front. 

For YSO RCW82 7, Fig. [B] shows again that H 2 is emitted 
over a wider region than the A"-band continuum. The H2 emis- 
sion is not centered on the YSO, but is offset to the west. In 
addition, diffuse emission concentrations are clearly visible es- 
pecially north-east of the YSO. The H2 velocity map constructed 
by fitting a Gaussian to the 2.12 fim line reveals that this blob 
has a radial velocity of -10 km s compared to -60 km s for 
the YSO Q. This may very well be the signature of a jet cen- 
tered on the YSO and interacting with molecular gas at the po- 
sition of the blob. At the distance of RCW 82, such a jet would 
have a projected length (YSO-blob) of ~0.06 pc. Besides this 
peculiar H2 emission, there is an arc-like H2 structure extending 
North and South of the YSO. Its velocity is comparable to that 
of the YSO and to that of the H 11 region (see above). It may thus 
correspond to the PDR on top of which the YSO is located. A 
similar morphology is identified in RCW120 49 (Fig. IB.5b and 
RCW120 103 (Fig. IB.6t . albeit with a lower significance (see 
Sect. IB.3) . Another convincing evidence for the presence of a 
jet comes from YSO RCW120 C4-67. It features a spectacular 
H2 map with an arc structure a few arcseconds away from the 
YSO (see Fig.fTBTl. Its radial velocity is comparable to that of the 
YSO, but the fact that the arc curvature is directed towards the 



4 The velocity map also reveal a region of lower velocity south-west 
of the YSO, but the emission is weaker here and the result is less sig- 
nificant. 
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YSO suggests that it might result from the interaction of a jet 
and molecular material. It may also be material ejected by the 
YSO in the past. 

From the analysis of the morphology and kinematics of a 
few sources, we conclude that, in some cases, we have been able 
to firmly establish the association of the YSO with the ioniza- 
tion front, strengthening the case of triggered star formation. 
We have also discovered the presence of a possible outflow in 
RCW82 7 traced by a H2 velocity gradient (and possibly also in 
RCW120 49 and RCW120 103). The geometry of the H 2 emis- 
sion in RCW120 C4-67 is also suggestive of the presence of a 
jet. 
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Fig. 14. Morphology of RCW120 51. Left panel: Bry emission 
map. Right panel: H2 emission map. The contours of the full 
K-band emission are overplotted. All images have been cre- 
ated from the initial datacubes smoothed with a spatial Gaussian 
(FWHM = 3 pixels). North is up, East is to the left. The axis 
indicate the offsets (in arscseconds) relative to the main YSO 
position (coordinates given in Tabled. 



6. Discussion 

6. 1 . Winds and bubbles 

The evolution of an H n region is governed by the ionizing 
radiation emitted by it s ionizing stars - s ee for example the 
first order analysis of iDvson &W illiams (1 19971) of the evo- 
lution of an H 11 region in an homogeneous medium. But 
it may also be influenced and even dominated by the action 
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Fig. 15. Morphology and dynamics of RCW82 7. Left panel: H2 
emission map together with contours of the full K-band emis- 
sion. Right panel: H2 2.12 pm velocity map (scale in km s~'). 
North is up, East is to the left. The axis indicate the offsets (in 
arscseconds) relative to the YSO position (coordinates given in 
Table©. 
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Fig. 16. Morphology of RCW120 C4-67. The image shows the 
H2 emission map (created from the initial datacubes smoothed 
with a spatial Gaussian - FWHM = 3 pixels). Contours of the 
full K-band emission are overplotted. North is up, East is to the 
left. The axis indicate the offsets (in arscseconds) relative to the 
YSO position (coordinates given in Table©. 



of the wind emitted by the central massive star dFrever et all 
12003b . It is often claimed that the many bubbles observed 
at all scales in the Gal actic plane result from the action 
of such powerful winds (Churc hwell et alj 120061) . While this 
is most likely the case when Wolf-Rayet stars are involve d 
dGarcia-Segura & Mac Lowlll995l:lBrighenti & D , Ercolelll997l) . 
the dominance of wind effects in the early evolution of interstel- 
lar bubbles remains to be established. This is even more true 
since the recent downward revision of the mass loss rates of 
O sta rs due to clumping (e.g. iBouret et alJfeOOSl: iFullerton et alj 

f200(k 

ICapriotti & Kozminskil (1200 ll) and iFrever et alj (120031) have 
developed analytical models and hydrodynamical simulations of 
H 11 regions with winds. At the beginning of its evolution, a 
young H 11 region should not be too much affected by the wind 
emitted by its ionizing star, if any. With time, the stellar wind 
gets stronger and eventually leads to the formation of a central 
cavity filled with very hot (T~ 10 6 K) low density, shocked gas. 
This gas emits at X-ray wavelengths. Such an extended X-ray 
emission is a strong indirect evidence for winds. It has been ob- 
served only in the dir e ction of three H n regions though: M17 
dTownslev et al.1 120031; iPovich et al.1 120071) . the Rosette nebul a 
(iTownslev et al.l2003l) . and the Orion nebula dGiidel etal.12 008). 
All regions lie relatively nearby (0.5 to 2 kpc), which might ex- 
plain that their X-ray emission is more easily detected than in 
distant regions. They are excited by a cluster of several O stars 
and have ages ranging from 1 (Orion) to about 3—4 Myr (Ml 7 
and Rosette). All three clusters contain early O stars (spectral 
types 04-06) with moderately strong winds (compared to early 
O type supergiants which have larger mass loss rates). This is 
similar to the three H 11 regions we have studied. We do not know 
if the extended X-ray emission is a common phenomena among 
bubbles or if it is due to the very massive stellar content of these 
regions. Unfortunately, RCW 79, RCW 82 and RCW 120 have 
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not been observed at X-ray wavelengths so the effects of winds 
in the H n regions cannot be quantified directly. 

Dust emission at 24 yum may be another indicator of the ef- 
fects of winds. Very small grains, heated by the stell ar radiation 
and o ut of thermal equilibrium, radiate at 24 //m dJones et alj 
1 19991) . As long as such emission is observed, it means that dust 
has not (yet) been swept-up by the stellar winds. The pres- 
ence of a central hole in Spitzer-MIPS 24 /mi emission maps 
of some H n regions might be an indication that winds are 
at work dWatson et alj 12008 *). Such a cavity could also be at- 
tributed to the radiation pressure of the central exciti ng stars or 
to dust destruction by the st ars' ionizing radiation dlnouell2002t 
iKrumholz & Matzneri 120091) . However Krumholtz & Matzner 
show that radiation pressure is generally unimportant for H n 
regions driven by a handful of massive stars (thus for our three 
regions). In any case if a shell is observed at 24 fim, it means 
that stellar winds have not removed all the dust from the H n re- 
gion. Fig.[l7]shows a composite image for the three H n regions 
studied here, with red, green and blue being respectively 24 fim, 
3.6/mi and Ha emission. The 24 fim emission traces the hot dust. 
In all three regions, one clearly sees the presence of 24 fim emis- 
sion bows surrounding the ionizing stars. The hot dust very close 
to those stars is thus either blown away or destroyed by the in- 
tense UV radiation. But 5" to 30" away, i.e. well within the bor- 
ders of the H ii regions, dust is detected. Hence, stellar winds 
have not (yet) removed all dust from the cavities, indicating that 
their effect on the dynamics of the H n regions is rather limited. 

To better assert the role of stellar winds and radiation in our 
objects, we focus on RCW 79. Since it hosts the largest number 
of ionizing stars of the three regions, it is the most appropriate 
to analyze the influence of stellar winds. 

The total number of ionizing photons emitted by the ionizing 
stars is 10 4982 s _I (see Table 14. 2T2l . This corresponds to an ion- 
izing luminosity of 10 39 31 erg s _1 . For an age of the population 
of 2.3 Myr, the ionizing energy released in the H n region is thus 
jq53.17 er g j n jjjg p rev i ous estimate we have assumed that the 
ionizing luminosity was constant with time. Since the stars are 
rather young and close to the main sequence, this is a reasonable 
approximation. According to the Geneva tracks, a 40 M star 
has T eff = 42560 K and log ^ = 5.34 on the ZAMS, and T eff = 

39170 K and log y- = 5.45 after 2.3 Myr. For the corresponding 
radius and using the number of ionizing photons per unit area 
(qo) from Ma rtins et alj (l2005al) . the values of Qo (total ionizing 
photon flux) are 10 4912 and 10 49 16 at and 2.3 Myr respectively. 
Hence our assumption is justified. 

The total wind mechanical luminosity at 2.3 Myr is 10 3618 
erg s _1 (for 12 O stars with a mass loss rate of 10~ 7 M Q yr~' 
and a terminal velocity of 2000 km s _1 ). Assuming again that 
this value is constant between and 2.3 Myr, one estimate a to- 
tal wind mechanical energy release of 10 50 04 erg. To check the 
validity of our as sumption, one can use the scaling relations of 
IVink et al.l d2001l) (their Eq. 24). Using the properties of a 40 M 
star on the ZAMS and after 2.3 Myr of ev olution (using once 
more the tracks of iMevn et & Ma ederll2005l) . one finds a mass 
loss rate difference of 0.3 dex. Hence, our wind mechanical en- 
ergy determination is at the very most overestimated by a factor 
2. 

The 2D simulations of lFrever et al] d2003l l2006h tackled the 
question of the effects of winds on the dynamics of H n regions. 
They showed that the presence of the wind triggered the forma- 
tion of structures ("fingers") in the H n regions, modifying the 
efficiencies of energy transfer (from stellar ionizing and wind 
mechanical energies to ionizing, kinetic and thermal energies in 



the H ii region). In their simulations, the ratio of ionizing to wind 
mechanical luminosities was of the order 100. They showed that 
in spite of this, winds could significantly affect the evolution of 
the H ii region. In our case, this ratio is at least 1000 (recall that 
our estimates of mass loss rates are upper limits). The main rea- 
son is that although the ionizing luminosities we find are similar 
to the ones used by Freyer et al., our mass loss rates are lower 
by nearly a factor of 10. The reason is the inclusion of clumping 
in our models. It is well known to lead to a dow nward revision 
of mass loss rates of O stars (e.g. lHillie^fl9 91). The values of 
M used by Freyer et al. were derived with homogeneous models 
and are thus higher than our determinations. With this increase 
by at least a factor of 10 in the ratio of ionizing to wind me- 
chanical luminosity, one might thus wonder whether winds still 
play a role. The evidence brought by the presence of 24 fim dust 
emission argues against a strong effect of stellar winds (at least 
in the first 2-3 Myrs). New dedicated simulations with reduced 
mass loss rates are encouraged to shed more light on this issue. 

6.2. YSOs: a possible spectroscopic evolution 

Among the YSOs observed on the borders of the three H n re- 
gions, a few have a TT-band spectrum entirely dominated by H2 
emission. We have seen in Sect. |5.3| that thermal emission due to 
shocks was favored to explain the different line ratios. Besides, 
most of these sources tend to be brighter and redder in the [24] 
vs A"-[24] diagram, indicating a larger amount of circumstel- 
lar material (Sect. 15.2b than for the other observed YSOs. All 
in all, this tends to favour a picture according to which the H2 
sources are in a relatively early phase of evolution. The large 
24 fim emission combined to the absence of direct indicator of 
the presence of a disk suggests that these sources still possess 
large envelopes. The H2 emission probably comes from either 
a jet-like structure interacting with this circumstellar material, 
or from an expanding shock created by stellar/disk wind and/or 
radiation pressure inside the envelope. One can further specu- 
late that objects showing both H2 lines and other features (Bry 
or CO bandheads emission) represent later stages of evolution in 
which the envelope partially disappeared, revealing the disk fea- 
tures. Their somewhat lower reddening and brightness in the [24] 
vs /T-[24] diagram lend support to this hypothesis. In that case, 
both Bry and CO emission could be produced in disks irradiated 
by the nascent star, Bry being emitted on the surface and CO 
close to the disk midplane, in zones shielded from UV flux. This 
could explain that the H2 lines are usually narrower than Bry in 
YSOs when both types of lines are detected (Table 0. The lat- 
ter being emitted on the disk's surface, it would be rotationally 
broadened, while the former would reflect the shock conditions 
in the dissolving envelope. We noted in Sect. l5.4l that a couple of 
objects showed radial velocities different for Bry and H2. In the 
suggested scenario, this could be explained if the H2 emission 
was produced by interaction of a jet and the envelope. If this pu- 
tative jet was directed towards us, one could see a blueshifted H2 
emission compared to the Bry line produced at the surface of the 
disk. 

Another piece of information comes from the probable de- 
tection of jets in a few sources (Sect. 15.41 ). All of these sources 
(RCW82 7, RCW 120 49, and RCW 120 C4-67) have been clas- 
sified as CO emitters. It means that they possess circumstellar 
disks. Since bip olar outflows are predicted to develop in accret- 
ing objects (e.g. IShu et al.ll 1 99*41) the detection of jets in several 
of our YSO confirms that we are witnessing accretion phenom- 
ena. Hence, those objects fit in the scenario presented above and 
correspond to evolutionary stages where the envelope has disap- 
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peared. This is consistent with the fact that all these sources are 
at the limit between class I and class II sources (see Fig. [10] and 
Table©. 

Obviously, the scenario presented above remains highly 
speculative. Future high spectral resolution spectroscopy of ob- 
jects showing CO, Bry and H2 emission will certainly provide 
more constraints on the geometry and dynamics of the line emit- 
ting regions. Nonetheless, if YSOs evolve spectroscopically ac- 
cording to the previous path and if they were formed simulta- 
neously, our observations indicate that in the three H 11 regions 
studied here YSOs evolve at different speed (presumably due to 
their different masses) since several types of objects are detected. 
Alternatively, if one drops the assumption that all YSOs formed 
at the same time, it means that the second-generation star for- 
mation event lasted less than 2 Myr. Said differently, the typical 
timescale for the evolution of YSOs (class to class II) should 
be < 2 Myr. 

7. Conclusion 

We have obtained near-infrared integral field spectroscopy of 
three H 11 regions (RCW 79, RCW 82, RCW 120) with SINFONI 
on the VLT. The main results of our study are: 

• We have identified the ionizing sources of all three H 11 re- 
gions. RCW 79 is powered by a cluster of O stars unraveled 
here for the first time. The most massive stars have a spectral 
type 04-6V/III. A number of later type OB stars are also 
detected. RCW 82 is ionized by two 09-B2V/III stars and 
a third star, probably of spectral type Be. Finally, a single 
O6-8V/III star ionizes the H 11 region in RCW 120. 

• The ionizing stars of each region have been analyzed with 
atmosphere models computed with the code CMFGEN. The 
derived stellar properties have been used to build the HR di- 
agram. In the case of RCW 79, the large number of stars 
allowed a reliable age determination: the ionizing cluster 
formed 2.3±0.5 Myr ago. For the other two H 11 regions, sim- 
ilar ages are suggested but are less well constrained due to 
the small number of sources. An upper limit of 10~ 7 M yr~' 
was derived for the mass loss rate of the ionizing stars of all 
three regions. 

• The cumulative mass loss rate due to the O stars at the cen- 
ter of RCW 79 is < 10~ 6 M yr -1 . The resulting wind me- 
chanical luminosity is about 10~ 3 the ionizing luminosity. 
This is ten times we aker than in the numerical simulations of 
Frev er et alj (|2003), raising the question of the quantitative 
role of winds on the dynamics of the three H 11 regions stud- 
ied here. The presence of hot dust emission close to the ion- 
izing stars argues against a strong influence of stellar winds. 

• Spectroscopy of the infrared sources on the borders of the 
H 11 regions revealed typical signatures of young stellar ob- 
jects. Four main categories of YSOs have been identified: (1) 
sources dominated by H2 emission lines, (2) sources show- 
ing mainly Bry emission, (3) sources with CO bandheads in 
emission and (4) sources with CO bandheads in absorption. 

• Near and mid-infrared color-color and color-magnitude di- 
agrams indicate that all YSOs have SEDs dominated by 
non stellar emission. In the Spitzer/IRAC [3. 6] -[4. 5] vs 
[5.8]-[8.0] diagram, essentially all YSOs are grouped at the 
position between class I and class II objects. Only a couple 
of H2 sources seem to be redder. The trend is more clearly 
seen in the [24] vs Zf-[24] diagram where the H2 sources 
form clearly a distinct group with larger 24 /im emission and 
redder /T-[24] colors. 



• The detailed analysis of the H2 emitting sources through line 
ratios and excitation diagrams indicate that thermal emission 
is the dominant mechanism in the majority of sources. Slow 
J-shocks or fast C-shocks are the preferred explanation for 
this emission. Some sources show a possible weak contri- 
bution of fluorescence emission, and three are dominated by 
this mechanism. 

• Measurement of radial velocities and line width was per- 
formed in a few sources. In general, the derived radial ve- 
locities are consistent with that of the ionized gas in the H 11 
region, confirming that the YSOs are well associated with the 
regions. The exception is RCW 120, for which most sources 
move 20 to 50 km s" 1 faster than the ionized gas. In sources 
showing both Bry and H2 lines, the former are systematically 
wider (150/250 km s~'vs non-resolved/100 km s~'), possi- 
bly due to a different spatial origin. 

• Morphological studies reveal that when both Bry and H2 
emission are present, Bry is centered on the position of the 
YSO while H2 emission is more extended. In some cases, 
the H2 emission peak is offset compared to that of the YSO. 
RCW 120 5 1 presents a double H2/Bry cone structure follow- 
ing the ionization front. RCW82 7 (and possibly RCW 120 
49 and RCW 120 C4-67) has a H 2 velocity structure typical 
of an outflow. 

• The spectrophotometric and kinematic properties of YSOs 
paint the following plausible picture: H2 dominated sources 
are stellar objects surrounded by an envelope in which 
shocks (due to jets or stellar/disk winds) produce H2 emis- 
sion; other YSOs correspond to more evolved sources in 
which the envelope partially or totally disappeared, reveal- 
ing a disk structure from which CO and/or Bry emission is 
produced. 

Future tailored observations of a few key YSOs will provide 
more information on the dynamics of the circumstellar material. 
Additional UV and optical spectroscopy of the ionizing stars of 
RCW 79 will help to refine the age estimate and to reveal a possi- 
ble age difference with the stars in the CHII region on the border 
of RCW 79. This would be a first direct proof of the existence 
of triggered massive star formation. Better estimates of the mass 
loss rates will also be possible, helping to quantify the role of 
stellar winds in the dynamics of H 11 regions. 
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Appendix A: Spectra 

In this Section, we present the spectra of our sources and the best 
fit to the ionizing sources (Figs. lATI to lA.4b . 

Appendix B: Notes on individual objects 

Here, we briefly summarize the main characteristics of each re- 
gion. Figs. IB.ll to IB. 71 provide source identification as well as 
linemaps of a few YSOs. 

B.1. RCW 79 

Eight sources were observed in this region. All present spec- 
troscopic features typical of YSOs. Sources 237, 91 and 3C3 
are resolved in multiple components (respectively Figs. lB.ll and 
IB~2l) . 

Sources 237-b and 1C4 show a decreasing continuum to- 
gether with Nai, Cai and CO absorptions. Source 1C4 also 
shows Bry in absorption. 

Sources 2C3, 3C3 and 237 are dominated by H2 emission. 
According to their GLIMPSE colors, they are Class I sources, 
3C3 and 237 being the reddest (high [3.6] - [4.5] color). This 
agrees with 3C3 and 237 having a rising continuum whereas 
2C3 has a flat one. Fig. IB. 31 shows the H2 emission of RCW79 
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Fig. A.2. SINFONI TT-band spectra of central stars of H n regions. Left: Be and late type stars in RCW 79. Right: ionizing stars of 
RCW82andRCW 120. 





Fig. A.3. Best CMFGEN fits (red) of RCW 79 ionizing sources spectra (black). 



2C3. While the Zf-band features a point source, a more elon- 
gated structure is seen in the H2 map. The emission peak in the 
latter does not coincide with the YSO in the Zf-band: a large 
amount of molecular material is still present around the YSO. 
We could not construct the H2 velocity map due to the weakness 
of the signal in the emitting regions. The H2 emission does not 
come directly from the YSO, and this should be kept in mind 
when interpreting the mixed results of Sect. 15.31 regarding the 
nature of the emission. 

Th e spectrum of 237 is very sim ilar to the one of source 
nr56c (Martin-Hernandez et al. 2008). A bright H2 emission fil- 
ament is observed in the field of source #237 (Fig. IB. lb . Its ve- 
locity is that of the H 11 region, and this H2 emission filament 



corresponds to the 8 yum emission that originates from the nearby 
PDR. 

Sources 18 (Fig.|EB, 91 (Fig.|R2]), 1 86 and 4C4 sho w wide 
Bry emission. A maser was identified by ICaswelll d2004l) in the 
region around source 91, but its association with any of the near- 
IR sources is not obvious (see ZA06). Sources 91,18 and 4C4 
are Class I sources, having similar colors in the GLIMPSE color- 
color diagram. 

The field of the compact H 11 region shows stellar sources 
and extended Bry and H2 emission (Fig. IB. 41 . The Bry emis- 
sion (ionized gas) is confined within a region bounded by a H2 
filament. OB stars are observed as dark spots in the Bry map. 
Surprisingly, the two brightest ionizing stars of the compact H 11 
region lie in the direction of the H2 filament and not in the ion- 
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Fig. A.5. Normalized SINFONI YSO spectra. Panel a, b, c and d show respectively the Bry, H 2 , CO emission and CO absorption 
dominated spectra. 



ized cavity; this might be a projection effect. The extended Bry 
emission and H2 filament have similar velocities, also similar to 
that of the large central RCW 79 region. 

Extended Bry and H2 emission are observed in the field of 
source 1C4 (not shown). The Bry emission is probably that of the 
ionized gas in RCW 79. The H2 emission comes from the PDR. 
Both emission indicate velocities similar to that of the central 
H 11 region. 

Source 4C4 is a tight cluster composed of one bright ob- 
jects and 6 much fainter components. They present an increasing 



spectrum and a faint Bry emission line. No extended emission is 
observed in this field. 



B.2. RCW 82 

Five candidates YSOs were observed in the direction of 
RCW 82, all situated on the border of this H 11 region, except 
possibly 51. Sources 4, 7, 35, and 51 have been classified as 
Class I by PO09. Source 98 has not been classified, due to 
uncertain J and H measurements and to an absence of 8 jim 
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measurement. 

Source 4 is almost featureless, with a faint and wide Bry 
line, and a rising spectrum. This agrees with a Class I source. A 
fainter source is present in the field, situated south-west of 4; it 
shows no features, except an increasing continuum. 

Source 7 has a rising spectrum with numerous H2 lines and 
CO headbands in emission. A disk is probably present. An H2 
jet is observed (Fig. |T3T>, confirming the youth of this object. Its 
evolutionary stage, I or II, is uncertain. 

Source 35 shows all types of signatures (Bry and H2 lines, 
CO emission, and an increasing spectrum). Its velocity field is 
unusual, with Bry emission at the velocity of the H n region, but 
strongly approaching H2 material. This source has probably a 
dominant disk and thus is probably a Class II object; ejection 
is possibly present, but needs confirmation. An extended 
H2 emission is observed in the field (see Fig. IB. II) . at the ve- 
locity of the H 11 region. It probably comes from the nearby PDR. 

Source 51 shows a strong Bry emission line, an increasing 
spectrum, and faint H2 and CO emission. Its high emission in the 
24 fim band indicates that it is probably a rather massive YSO. 
It lies at the center of a cluster (Fig. IB. It : several faint con- 
tinuum sources are present in the field (see also Fig. 17 in PO09). 

Source 98 (Fig. IB.lt shows faint CO features in absorption; 
however it is most probably a young source (and not an evolved 
star) as it presents an increasing spectra. A bright extended 
H2 emission is present over the whole field (Fig. A. 8), at the 
velocity of the H 11 region. Most probably it is emitted in the 
vicinity of the PDR as confirmed by the line ratios indicative of 
a non-thermal process (^21 ~ 3.5 and „ 10) 

1 '2-15(1) '3-25(3) 



B.3. RCW 120 

Ten sources were observed in the direction of RCW 120. Eight 
of them present spectroscopic features characteristic of YSOs. 
Two of them are different. They are: 

• RCW120 39 was classified as a Class II in DE09. The 
SINFONI spectrum shows that it is most probably an evolved 
star, as it presents strong CO absorption features and no other 
signature typical of YSOs. 

• RCW 120 C4-64 was also classified as an intermediate 
Class I/II source in DE09. However its spectrum is slowly de- 
creasing; there is no Bry nor H2 emission lines; only the CO 
bands are observed in absorption. Thus it is probably also an 
evolved star, located by chance on the border of the C4 conden- 
sation. 

Note that RCW 120 lies only -12° of the Galactic center, 
and thus evolved stars of the Galactic bulge may be present in 
the background. 

RCW120 76 (Fig. ED, 106 and 107 are rather similar, 
almost featureless, with no or faint Bry line, no or faint H2 lines 
and no CO emission bands, along with an increasing continuum. 
They are possibly the less evolved of the YSOs surrounding 
RCW 120, mostly dominated by their envelope. They were 
classified as respectively Class I-II, Class I, and Class I-II in 
DE09, in good agreement with the present results. 



RCW 120 49 and 103 are similar, with strong Bry emis- 
sion, numerous H2 lines and CO bands, superimposed on an 
increasing continuum. They are possibly more evolved than the 
preceding group, the presence of strong CO emission bands 
pointing to a dominant disk. They were classified as Class I-II 
and Class I in DE09. 

Fig. IB. 51 shows the H2 emission and velocity maps of 
RCW120 49. The H2 emission shows multiple peaks none of 
which is located at the A'-band position of the YSO. The H2 
velocity map reveals that the two H2 peaks on each side of the 
YSO have different radial velocities. We note however that the 
H2 material at velocities of about -80 km s _1 corresponds to a 
weak emission region (it is slightly offset compared to the emis- 
sion peak seen in the H2 map). The H2 profile at this position 
is also peculiar in the sense that it probably results from multi- 
ple components, being rather broad and featuring a "shoulder" 
on the blue side of the main emission peak. On the other hand, 
the bright emission peak south-east of the YSO is significantly 
detected at a lower velocity than the YSO itself and do not show 
any peculiar line profile. The analysis of the H2 emission lines in 
Sect. |5.3| indicated thermal emission with a possible small contri- 
bution from non thermal emission^. The presence of a jet would 
explain both this emission and the velocity map. 

We have seen previously that YSO RCW 120 103 had pe- 
culiar kinematical properties, with H2 moving much faster than 
Bry (Sect. l5~4l . Fig. lB.6l shed more light on this problem. The H2 
emission peaks at the same position as the /f-band emission. But 
an arc/spiral is detected south of the YSO, showing that molecu- 
lar gas is present on wide scales. The velocity map confirms that 
a minimum (around -90 km s~') is located on the YSO. But this 
map tend to show that a velocity gradient exists on a direction 
south-east/north-west across the position of the YSO. A similar 
map in Bry did not reveal any structure. We note however that 
the velocity difference between both sides of the H2 structure is 
of the order of 20 km s _1 , just at the limit of the resolution of our 
velocity measurements. Hence, we refrain from interpreting too 
much this possible structure and simply conclude that the geom- 
etry and dynamics of RCW 120 103 is probably rather complex. 
This is confirmed by the velocity of the extended H2 emission 
south of the YSO, which is intermediate between the velocity of 
the H 11 region and that of the YSO. We also note that RCW82 
35 which showed the same velocity discrepancy as RCW120 103 
(H2 velocity much larger than Bry velocity) did not reveal any 
velocity structure. 

The field of RCW120 51 contains several sources (Fig. [T4b. 
all observed in the direction of the ionized gas. We see on the 
spectra the Bry emission of the ionized gas. The three brightest 
sources have a decreasing continuum; source 51b and 51c 
present CO absorption bands. These two sources are probably 
stars and not YSOs. Source 51a is probably an YSO, as its 
spectrum show Bry emission and no CO absorption. However 
we do not confirm the Class I classification of DE09. YSO 
51a lies at the tip of a "finger", a structure of the ionization 
front bright in the mid-IR. The proximity of this structure 
has probably hampered the photometric measurements of this 
object. 

Condensation 4 contains several sources with different 
properties (Fig. IB. 7b . Sources A and Ba/Bb are respectively 
at the center of zones of extended 8 //m and 24 fim emission 



5 Note that the analysis was performed on the spectrum extracted at 
the position of the IT emission peak. 
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resembling small PDRs. The JHK 2MASS photometry of the 
central objects point to late B stars (DE09). The SINFONI 
spectra of these 3 sources show decreasing continua. Source 
A shows both Bry and CO absorption. This might result from 
the blending of two sources, one with each type of features, 
or be indicative of a medium type star. Source Ba has Bry in 
absorption and is most likely a B star, while source Bb features 
CO absorption bands and is either a late type star or a YSO (see 
Sect. 15. U . Three sources are observed in the field of RCW120 
69. They all have a flat continuum. RCW120 C4-69a is fea- 
tureless. RCW120 C4-69c has CO absorption bands. RCW120 
C4-69b show faint CO absorption bands but Bry in emission. 
Extended H2 emission is observed south-east of the main 
source. These sources are probably evolved YSOs, RCW120 
C4-69b being the youngest. RCW120 C4-69, considered as a 
single source, was classified Class I by DE09. However it is not 
a Class I, but more probably a Class II or III. 




Fig. 17. Composite color image around the ionizing stars or 
cluster in the three H 11 regions: red is the Spitzer-MIPSGAL 
emission at 24 fim, green is the Spitzer-GLIMPSE emission at 
3.6 jjm, and blue is the Ha emission. The black dashed line en- 
closes the zones where the 24 fim emission is saturated. North is 
up and east is left. The ionizing star(s) or cluster lie in the centre 
of the H n regions, in a zone void of dust emitting at 24 fim. They 
are indicated by blue arrows. 
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Fig.A.l. SINFONI /f-band spectra of the ionizing sources of 
RCW 79: O stars. The main stellar features are indicated. 
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Fig. A.6. Spectra of YSOs with composite features. 




Fig.B.3. Morphology of RCW79 2C3. The image shows the H 2 
emission (created from the initial datacubes smoothed with a 
spatial Gaussian - FWHM = 3 pixels). North is up, East is to 
the left. The axis indicate the offsets (in arscseconds) relative to 
the YSO position (coordinates given in Tabled. 
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Fig. A.4. Best CMFGEN fits (red) of RCW 82 and RCW 120 
ionizing sources spectra (black). 
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Fig. B.4. Bry (top) and H2 maps of the the compact H n region 
on the border of RCW 79. The contours of the full K-band emis- 
sion are shown in black. North is up and East is to the left. The 
axis indicates the RA and DEC offsets (in arcseconds) relative to 
source CHII #6. 



Fig.B.5. Morphology and dynamics of RCW120 49. Top: H 2 
emission map with contours of the full K-band emission over- 
plotted. Bottom: H2 velocity map (scale is in km s -1 ). North is 
up, East is to the left. The axis indicate the offsets (in arscsec- 
onds) relative to the YSO position (coordinates given in Table 
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Fig. B.7. Identification of YSOs in condensation 4 of RCW 120. 
North is up and East is to the left. The axis indicates the RA and 
DEC offsets (in arcseconds) relative to source 67. 
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Fig. B.l. Top panels: H2 maps of YSOs with contours of the full K-band images overplotted. The H2 emission is more extended 
than the bulk of the K-band emission. Bottom panels: Bry and H2 maps of RCW82 35 and RCW79 18 together with contours of the 
full K-band emission. In all panels, North is up and East is to the left. The axis indicates the RA and DEC offsets (in arcseconds) 
relative to the main YSO (see Table[3]for coordinates). 
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Fig.B.2. Identification of the sub-components of YSOs RCW79 91 (top panels) and RCW79 3C3 (bottom panels) together with 
Bry and H2 maps. North is up and East is to the left. The axis indicates the RA and DEC offsets (in arcseconds) relative to the main 
YSO (see Table[3]for coordinates). 



